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Preface 

This  thesis  is  a  continuation  and  extension  of  previous 
work  by  graduate  students  at  the  Air  Force  Institute  of 
Technology  in  the  area  of  parameter  estimation  using  order 
statistics  of  samples  from  a  given  distribution.  The  tables 
of  linear  coefficients  developed  in  this  report  will  enable 
the  user  to  obtain  the  best  linear  invariant  estimate  of 
the  location  and  scale  parameters  of  the  Cauchy  distribution 
for  sample  sizes  of  N=5(l)20  very  efficiently.  An  attempt 
was  made  in  the  report  to  provide  a  clear  development  of  the 
theory  by  which  these  linear  coefficients  are  obtained.  In 
addition,  the  Fortran  program  required  to  calculate  and  table 
the  linear  coefficients  is  included  in  Appendix  C.  The  sub¬ 
routine  used  to  solve  the  matrix  equations  is  a  modification 
of  the  Matrix  Equation  Solver  Fortran  Subroutine  from  the 
Computer  Science  Center,  Wright-Pattcrson  AFB,  Ohio. 

I  also  wish  to  acknowledge  my  debt  to  Professor  Albert 
fl.  Moore,  my  advisor,  for  proposing  this  area  of  study  and 
for  his  assistance  and  encouragement. 


Ralph  M.  Spory,  Jr. 
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Linear  coefficients  which  can  be  applied  to  sample  data  from,  a  Cauchy 
distribution  to  obtain  estimates  of  tho  location  and  scalo  parameters  are 
developed  and  tabled.  Several  previous  works  have  presented  such  tables  for 
nearly  best  linoar  unbiased  estimation  and  beet  linear  unbiased  estimation  of 
tho  parameters.  The  estimates  developed  in  thia  paper  aro  best  in  tho  sense  that 
they  possess  minimum  moan  scuaro  error.  By  using  oxact  values  of  tho  moans, 
variances,  and  covariances  of  the  Cauchy  standardised  order  statistics  and  min¬ 
imising  tho  mean  square  error  function,  matrix  equations  aro  developed  and 
solved  to  obtain  the  required  coefficients.  These  coefficients  and  values  of 
the  MSB  are  tabled  for  minimally  consorecl  sample  sines  of  5  to  20  and  for  sam¬ 
ples  which  have  been  additionally  censored  from  abovo  and  symmetrically. 
Procedures  for  using  the  tobies  and  several  illustrative  calculations  demonstrate 
the  simplicity  of  this  estimation  technique.  Tho  Fortran  programs  required  to 
calculate  and  tablo  the  above  values  aro  included  in  Appendix  C. 
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Abstract 

Linear  coefficients  which  can  be  applied  to  sample  data 
from  a  Cauchy  distribution  to  obtain  estimates  of  the  loca¬ 
tion  and  scale  parameters  are  developed  and  tabled.  Several 
previous  works  have  presented  such  tables  for  nearly  best 
linear  unbiased  estimation  and  best  linear  unbiased  estima¬ 
tion  of  the  parameters.  The  estimates  developed  in  this 
paper  are  best  in  the  sense  that  they  possess  minimum  mean 
square  error.  By  using  exact  values  of  the  means,  variances, 
and  covariances  of  the  Cauchy  standardized  order  statistics 
and  minimizing  the  mean  square  error  function,  matrix  equa¬ 
tions  are  developed  and  solved  to  obtain  the  required  coef¬ 
ficients.  These  coefficients  and  values  of  the  MSE  are 
tabled  for  minimally  censored  sample  sizes  of  5  to  20  and 
for  samples  which  have  been  additionally  censored  from  above 
aud  symmetrically.  Procedures  for  using  the  tables  and 
several  illustrative  calculations  demonstrate  the  simplicity 
of  this  estimation  technique.  The  Fortran  programs  required 
to  calculate  and  table  the  above  values  are  included  in 
Appendix  C. 
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CONDITIONAL  BEST  LINEAR  INVARIANT  ESTIMATION  OF  THE 
LOCATION  AND  SCALE  PARAMETERS  OF  THE  CAUCHY 
DISTRIBUTION  BY  THE  USE  OF  ORDER  STATISTICS 


I.  Introduction 


Statement  of  the  Problem 

Objective.  The  objective  of  this  thesis  is  to  develop 
a  table  of  linear  coefficients  which  can  be  easily  applied 
to  sample  data  from  a  Cauchy  distribution  to  determine  esti¬ 
mates  of  the  location  and  scale  parameters  of  the  distribu¬ 
tion.  These  estimates  of  the  parameters  are  conditional 
best  linear  invariant  estimates.  These  terms  and  the  proper¬ 
ties  of  the  estimators  are  defined  in  the  next  section. 

Definition  of  Terms.  The  Cauchy  distribution  is  a  con¬ 
tinuous  distribution  which  is  frequently  introduced  to  stu¬ 
dents  as  an  example  of  a  distribution  for  which  the  moments 
do  not  exist  (Ref  6:134).  The  cumulative  distribution 
function  (cdf)  is  given  by 


F(x)  =  1+1  ARCTAN  ^  ,  -»<x«*  m 

2  it  s  *  v  ' 

and  the  probability  density  function  (pdf)  is  given  by 


f(x) 


_ s _ 

Tf[s2+(x-t)2] 


-cu<x<t» 


(2) 


where  s  is  the  scale  parameter  and  t  is  the  location  param¬ 
eter. 
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Conditional  estimation  of  a  parameter  is  defined  as 
estimation  of  an  unknown  parameter  when  the  second  parameter 
is  known.  In  this  case  the  location  parameter  is  estimated 
conditioned  on  the  value  of  the  scale  parameter,  and  the 
scale  parameter  is  estimated  conditioned  on  the  value  of 
the  location  parameter. 

The  best  linear  invariant  estimator  is  best  in  the 
sense  that  it  is  a  minimum  mcan-square-error  estimator.  The 
mean  square  error  of  the  estimator  is  given  by 

MSE  =  E[(s*-s)2]  (3) 

where  s  is  the  true  value  of  the  parameter  and  s*  is  the 
estimated  value  of  the  parameter. 

For  conditional  estimation  the  mean  square  error  is 
given  by 

MSE  =  E[(s*|t-s)2]  (4) 

where  s*|t  is  the  estimator  of  the  parameter  conditioned  on 
the  value  of  t,  and  t  is  the  value  of  the  second  parameter. 

Linear  estimation  is  a  technique  based  on  the  theory  of 
order  statistics,  where  each  ordered  sample  value  is  assignee! 
a  weight,  or  linear  coefficient,  and  the  coefficients  are 
calculated  so  as  to  obtain  the  best  estimate  of  the  parameter. 
In  this  case  the  best  estimate  is  the  invariant  estimate  as 
defined  above.  If  one  uses  this  method,  the  estimate  of  the 
parameter  is  given  by 
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n 

s*  =  E  ajXi  (5) 

i=l 

where  the  a-['s  are  the  linear  coefficients,  the  Xj/s  are 
the  ordered  sample  values  or  the  order  statistics,  and  n  is 
the  sample  size. 

The  conditional  estimator  is  then  given  by 

n 

s*|t  **  E  aiXi-At  (6) 

i=l 

where  t  is  the  known  parameter  and  A  is  a  constant  that  is 
determined  by  the  calculation.  The  significance  of  this 
constant  will  be  explained  in  Chapter  IV. 

Order  statistic  theory,  the  Cauchy  distribution,  and 
the  method  used  to  solve  for  the  linear  coefficients  are 
developed  further  in  the  remaining  chapters  of  this  paper. 

Significance .  When  the  location  and  scale  parameters 
of  the  Cauchy  distribution  are  known,  the  function  is  com¬ 
pletely  defined,  and  it  can  be  used  in  a  decision-making 
process  where  one  is  working  with  data  which  is  distributed 
according  to  the  Cauchy  law.  The  tables  of  coefficients 
developed  ii.  this  thesis  will  allow  the  user  to  estimate  the 
values  of  these  unknown  parameters.  These  estimates  can  be 
obtained  very  efficiently  with  only  the  use  of  a  desk  cal¬ 
culator.  The  ordered  sample  values  arc  simply  multiplied 
by  the  appropriate  coefficients  and  the  results  summed  to 
calculate  the  estimate. 
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This  method  of  estimation  provides  a  great  savings  in 
time  in  the  case  of  the  Cauchy  distribution,  as  the  tradi¬ 
tional  methods  of  obtaining  these  estimates  either  do  not 
apply  to  the  distribution  or  are  very  time-consuming.  Two 
of  these  methods  are  the  method  of  moments  (Ref  11:186)  and 
the  method  of  maximum  likelihood  (Ref  11:170).  The  method 
of  moments  cannot  be  applied,  since  the  moments  of  the 
Cauchy  distribution  do  not  exist,  and  the  maximum  likeli¬ 
hood  estimate,  which  is  convenient  to  obtain  for  some  other 
distributions,  requires  a  great  amount  of  computational 
effort.  Barnett  (Ref  1)  points  out  that  the  frequent  occur¬ 
rence  of  multiple  zeros  of  the  derivatives  of  the  logarithm 
of  the  likelihood  function  requires  a  complete  scan  of  the 
likelihood  function  to  locate  the  maximum  which  corresponds 
to  the  maximum  likelihood  estimate.  The  tables  of  linear 
coefficients  in  this  thesis  will  provide  the  user  with  an 
estimate  of  these  parameters  for  sample  sizes  of  5  to  20. 

Background  Information 

Work  on  parameter  estimation  based  on  order  statistic 
theory  has  been  carried  out  by  the  students  at  the  Air  Force 
Institute  of  Technology,  under  the  direction  of  Professor 
Albert  H.  Moore,  and  sponsorship  of  Dr.  II.  Leon  Harter  (ARL, 
Wright-Patterson  AFB) ,  since  1963.  These  works  include 
parameter  estimation  of  the  Cauchy,  Weibull,  normal,  log¬ 
normal,  logistic,  and  extreme  value  distributions.  Parameter 
estimation  of  the  Cauchy  distribution  includes  the  works  of 
Chamberlain  (Ref  2),  Jonson  (Ref  9)  and  Stark  (Ref  16). 
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Chamberlain  developed  and  tabled  the  coefficients  for  nearly 
best  linear  unbiased  estimation  of  the  location  and  scale 
parameters  for  sample  sizes  15(1)40.  Jonson  computed  the 
coefficients  for  conditional  best  linear  unbiased  estimation 
of  the  parameters  of  the  Cauchy  distribution  and  compared 
the  efficiency  of  these  estimators  with  the  efficiency  of 
Stark's  best  linear  unbiased  estimators.  The  estimators 
developed  by  Jonson  and  Chamberlain  are  called  nearly  best 
estimators  because  the  approximate  values  of  the  covariances 
of  the  order  statistics  given  by  Blom's  approximation  were 
used  instead  of  the  exact  covariances  of  the  order  statis¬ 
tics.  Stark's  work  developed  linear  coefficients  for  simul¬ 
taneous  estimation  of  the  location  and  scale  parameters  by 
using  the  exact  values  of  the  means,  variances,  and  covar¬ 
iances  of  the  standardized  order  statistics. 

The  works  of  Chamberlain  and  Jonson  are  based  largely 
on  the  methods  presented  by  Barnett  (Ref  1:1205).  Barnett 
tabled  the  coefficients  for  the  best  linear  unbiased  esti¬ 
mate  of  the  location  parameter  of  the  distribution  for 
sample  sizes  of  5  to  20.  The  exact  values  of  the  neans, 
variances,  and  covariances  of  the  Cauchy  order  statistics 
were  calculated  to  four  decimal-place  accuracy.  The  values 
of  the  covariances  were  obtained  by  numerical  integration 
of  expressions  containing  the  joint  pdf  of  the  order  sta¬ 
tistics  (see  Chapter  IV).  These  functions  were  integrated 
over  the  relevant  triangular  region  by  a  two-dimensional 
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extension  of  the  composite  Simpson  procedure.  Although  this 
computation  required  a  large  number  of  steps  to  obtain  the 
desired  accuracy,  it  did  prove  feasible. 

The  median  of  the  sample  data  has  traditionally  been 
used  as  an  estimate  of  the  location  parameter  of  the  Cauchy 
distribution.  Cramer  (Ref  4:708)  states  that  the  variance 
of  this  estimator  is  *2/4n  for  large  samples.  Rider  (Ref 
13:322)  shows  that  this  is  not  a  very  accurate  estimate  of 
the  variance  of  the  median  for  small  sample  sizes.  Rider 
has  tabled  the  actual  variances  of  the  median  for  small 
sample  sizes. 

In  1964,  Rothenberg  et  al.  proposed  a  class  of  esti¬ 
mators  of  the  location  parameter  of  the  Cauchy  distribution 
which  is  the  arithmetic  average  of  a  central  subset  of  the 
sample  order  statistics.  The  sample  median  is  a  member  of 
this  subset,  but  it  was  shown  that  the  average  of  the  middle 
quarter  of  the  ordered  samples  has  a  lower  asymptotic  var¬ 
iance  than  does  the  median. 

In  1970,  Chan  (Ref  3:851)  proposed  a  conditional  asymp¬ 
totically  best  linear  estimator  of  the  location  and  scale 
parameters  based  on  a  few  of  the  order  statistics.  He  has 
tabled  coefficients  for  K=l(l)10  where  K  is  the  number  of 
order  statistics  selected  from  a  large  sample.  These  esti¬ 
mators  yield  more  than  92  percent  asymptotic  relative  ef¬ 
ficiency,  in  the  Craraer-Rao  sense,  for  K>4. 
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Report  Orga^. ;  zation 

In  order  to  develop  the  linear  coefficients  for  the 
conditional  best  linear  invariant  estimation -of  the  param¬ 
eters  of  the  Cauchy  distribution,  the  distribution  and 
methods  by  which  it  is  generated  will  be  discussed  in 
Chapter  II.  Chapters  III  and  IV  present  order  statistic 
theory  and  the  linear  estimation  procedure.  Chapter  V 
describes  the  tables  of  linear  coefficients  and  gives  exam¬ 
ples  of  the  calculations  required  to  obtain  the  desired 
estimators.  Tables  of  the  values  of  the  mean-square-error 
function  and  the  linear  coefficients  are  included  in  Ap¬ 
pendices  A  and  B.  Appendix  C  contains  the  Fortran  programs 
required  to  calculate  and  table  the  above  values. 

Assumptions 

There  are  two  assumptions  made  in  this  report.  It  is 
assumed  that  the  sample  data  are  known  to  come  from  a 
Cauchy  distributed  parent  population  and  that  the  parameter 
not  being  estimated  is  knoim,  in  the  case  of  conditional 
estimation.  In  the  case  of  simultaneous  estimation  it  is 
only  assumed  that  the  sample  data  are  known  to  come  from  a 
Cauchy  distributed  parent  population. 
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II.  The  Cauchy  Distribution 

Introduction 

The  Cauchy  distribution  is  a  continuous,  symmetric  dis¬ 
tribution.  The  cdf  and  pdf  are  given  by  Equations  (1)  and 
(2).  The  plot  of  the  pdf  is  similar  to  that  of  the  more 
familiar  normal  distribution,  except  that  the  curve  ap¬ 
proaches  the  axis  much  more  slowly  and  the  tails  are  thicker. 
Figure  1,  on  the  following  page,  is  a  plot  of  the  density 
function  for  three  different  scale  parameters.  Feller 
(Ref  5:57)  provides  an  excellent  description  of  the  Cauchy 
distribution,  its  peculiarities,  and  methods  by  which  the 
distribution  is  generated. 

Example  of  the  Cauchy  Distribution 

A  mirror  is  arranged  parallel  to  an  opposing  wall  at  a 
distance  S  from  the  wall,  and  the  mirror  is  free  to  rotate 
on  a  vertical  axis  at  A  which  is  located  on  a  line  per¬ 
pendicular  to  the  wall  at  0.  The  angle  0  is  measured  from 
this  line  to  the  perpendicular  to  the  surface  of  the  mirror. 
The  mirror  reflects  a  ray  of  light  on  the  wall  at  a  distance 
X  from  the  point  0.  Now,  if  the  angle  0  is  chosen  at  random 
between  -ir/2  and  tt/2,  the  random  variable,  X,  is  Cauchy 
distributed  and  the  density  of  the  distribution  is  given 
by  Equation  (2)  with  a  location  parameter  of  zero  (Ref  5:57). 

i’he  density  of  the  random  variable,  X,  can  easily  be 
verified  by  a  method  due  to  Meyer  (Ref  10:88-89).  In  the 
above  example 
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f((0  -  I  for  -j<0<j 


is  a  uniform  density.  Then 


g(x) 


1  d0 
ir  dx 


(7) 


where  0  =  tan'1  (£) .  Therefore 


gW  -  i 

—K  1 
i+£:  s 

s^ 

(8) 

g«  -  i 
11 

[S2.x2]  for 

(9) 

This  is  the  Cauchy  pdf  with  scale  parameter  s  and  location 
parameter  zero. 


Generation  of  the  Cauchy  Density 

The  Cauchy  density  may  be  generated  in  many  ways.  The 
Student’s  t  density  with  n=l  is  a  Cauchy  density.  If  x  and 
y  are  two  independent  random  variables  from  a  standardized 
normal  distribution,  the  quotient  of  these  random  variables 
is  a  standardized  Cauchy  distribution.  In  addition,  if  the 
random  variable  x  is  Cauchy  distributed  with  a  scale  par¬ 
ameter  of  1  then  1/x  has  the  same  density.  Once  again, 
these  densities  may  be  verified  by  the  method  described  by 
Meyer  (Ref  5:109-110). 


Properties  of  the  Cauchy  Distribution 

Due  to  the  thick  tails  of  the  Cauchy  distribution,  esti¬ 
mation  of  its  center  is  very  difficult.  The  moments  of  the 
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Cauchy  distribution  de^not^exist.  Jonson  (Ref  9:11)  shows 
the  characteristic  function  of  the  Cauchy  random  variable 
to  be 

Cx(t)  s  exp (ict-b| t | )  (10) 

where  C  is  the  location  parameter  and  b  is  the  scale  par- 
aic^ter,  and  that  the  moments  do  not  exist  since  the  partial 
derivatives  of  Cx(t)/i^  with  respect  to  t  evaluated  at  t=0 
are  infinite. 

It  can  also  be  shown  directly,  that  the  first  moment 
about  the  origin  of  the  Cauchy  distribution  does  not  exist 
(Ref  6:145). 


11 


GAM/MATU/72-3 


III.  Order  Statistic  Theory 

Introduction 

To  calculate  the  conditional  best  linear  invariant 
estimates  of  the  parameters  of  the  Cauchy  distribution  the 
exact  values  of  the  means,  variances,  and  covariances  of 
the  Cauchy  order  statistics  are  required.  The  order  sta¬ 
tistic  theory  and  application  of  this  theory  to  the  Cauchy 
distribution  will  be  reviewed  in  this  chapter.  Reference  15 
is  a  rather  complete  collection  of  contribut ’cna  to  order 
statistic  theory  and  includes  articles  up  to  1962.  The 
following  development  follows  that  of  Chapter  II  from  the 
above  reference. 

Older  Statistics 

Definition.  If  a  random  sample  of  size  n  (xj,  X£,  ... 
Xn)  is  taken  from  a  population,  these  independent  random 
variables  can  be  rearranged  so  that 

X(l)iX(2)i*“lx(n) 

When  the  variables  arc  arranged  in  order  of  magnitude,  they 
arc  called  order  statistics  of  the  sample.  Since  these 
samples  are  from  a  continuous  population 

P(x(i)  »  x(j)'=!  0  for  a11  i^j 

Density.  The  pdf  of  the  ith  order  statistic,  x(i)>  is 
given  by 
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gCx(i))  =  (i-i) r(n-i) ! fpCx(i) J1  1[1_FCx(i)3n  lfCx(i)) 

(ID 

where  F(x^p  =  cd£  of  x  evaluated  at  x~x^ 

f(x(ip  c  pdf  of  x  evaluated  at  x=X(jj  (Ref  15:12) 

The  joint  distribution  of  the  ith  and  jth  order  statistics 
(i<j)  is  given  by 

8(x(i)  >x(j) )  =  ]jrrjT ( f)T (n -ITT ^ F (i ) 5 J 1  1fF(x(j)) 

(12) 

F(x(i)))^"1"1* [l~F(x(j))]n"jf(x(i))f(x(j)) 

for  x(i)<x(j)  (ReF  15:12) 

From  Equations  ^11)  and  (12)  expressions  for  the  expected 
values  and  covariances  of  the  Cauchy  order  statistic  can  be 
developed. 

Expected  Values.  Let  x  be  a  continuous  random  varia¬ 
ble  with  pdf  f,  then  the  expected  value  of  x  is  given  by 

00 

E(x3  «=  fx f(x)dx  (Ref  10:121)  (13) 

•00 

Using  Equations  (13)  and  (14) ,  one  finds  that  the  expected 
value  of  the  ith  order  statistic  is  given  by 

00 

Elx(i)3  85  /x(i)S(x(i))dx(i)  (14) 

-•0 

And  the  expected  value  of  the  product  of  the  ith  and  jth 
order  statistics  is  given  by 
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CO  CO 

E[x(1)xa)]  =  /  /x(i)xa)g(*(i)*0))dx(ijdx03  (15) 

—  CO  —  oo 

where  g(x^j  ,x^ j -j )  is  as  defined  by  Equation  (12). 

Variance  and  Covariance.  A  two-dimensional  random 
variable  Cx(i)*x(j))  possesses  a  property  called  the  covar¬ 
iance  of  x(i)>x(j)*  In  this  case  the  random  variables  are 
order  statistics,  and  in  a  sense  the  covariance  is  a  measure 
of  the  dependence  between  the  two  values  of  the  order  sta¬ 
tistics.  The  covariance  of  (x(i)»x(j))  is  formally  defined 
as  the  product  moment  about  the  respective  expected  values 
of  the  order  statistics.  Meyer  (Ref  10:144)  defines  the 
covariance  of  the  two  random  variables  as  follows: 

Cov(X(jj  >xQp  =  E{  [X(i)-E(X(i))]  [X(j)-E(-x(j))]}  (16) 

It  can  be  easily  shown  from  (16)  that: 

Cov(x(i)  ,x(j))  »  Etx(i)x(j)l“Efx(i)]  •Etx(j)]  (17) 

The  variance  of  the  ith  order  statistic  may  be  considered  as 
a  special  case  (i=j)  of  Equation  (16)  where  the  Var  lx(i)] 
is  given  by 

Var(x(i)]  *  E{  [X(i)-E(x(i)]2}  (18) 

With  the  expressions  presented  in  this  chapter,  the  expected 
valves,  variances,  and  covariances  of  the  Cauchy  order  sta¬ 
tistics  can  be  developed. 

Standardized,  Cauchy  Order  Statistics.  Let  x^jj  ,  x(2) > 
,,,x(n)  Ee  a  sct  ordered  sample  values  from  a  Cauchy 
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distributed  parent  population.  A  new  set  of  standardized 
order  statistics  can  be  developed  where 

-«<U(i)<~  (19) 

with  Uqj  <***<U(n) 

_ s _ 

If  the  pdf  of  x  is  f(x)  =  v [s2+(x_t) 2]  (2) 

then  the  density  of  U  is  given  by 

P(u)  =  -»<U<»  (20) 

and  its  cdf  is  given  by  P(u)  =  \  +  -  ARCTAN  u  (21) 

The  pdf  of  the  ith  standardized  order  statistic  from  Equa¬ 
tion  (11)  is  given  by 

-  (i-l)"(n-'iyrt‘,Cu(i)]1-1[l-P(u(i)]n-ip(u(1j)(22) 

and  the  joint  pdf  of  the  standardized  order  statistics  is 
given  by  Equation  (12)  with  x  replaced  by  u,  g  by  p,  F  by  P 
and  f  by  p. 

The  above  expression  can  be  simplified  by  making  the 
following  substitution: 

Let  =  ARCTAN  «(i) ,  -  "C0<u(i)<00  (23) 

where  du(jj  =  (1  +  tan20^)d0j  (24) 

and  =  ARCTAN  u(j)  -*<0<£  (25) 

then  from  Equation  (22) 
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q(UU)  *  Ci-l)l(n-i)!,n  [6  -  £] ^  [f-e]""1 (26) 


q(u(i)u(j)5  =  (i-1)  !  (j-i-1)!  (n-j)  ! 7in  [0+IJ 


(27) 

[(8+1)-  (0+i) ]5"i_1  •  ri-0in'j - ~ - ^ - 

%  2  2  (l+tan^Q)  (l+tan^#) 

By  substituting  these  equations  into  Equations  (14),  (IS) , 
(17),  and  (18)  the  desired  expressions  for  the  expected 
values,  variances  and  covariances  of  the  standardized  order 
statistics  are  obtained. 


E[“a)]  ’  (TTJUn-ijM  |tane[e+|]1'1(J-0]n"ide: 


2 

Cov[uci),u(j)]  -  (TJjjjy.JllJKn -j)ltTr  //tan0tane[e+’]i-1 

_TT  _ir 

7  7 

•  [  («♦$)  -  C0+i)  P  ’i'1  d9d0-- - - - - 

22  2  (i-l)!(n-i)!(j-l)!(n-j)Un 

IT  £ 

• /tan9  te+f]1'1 [f-9]n"id9/tan? [0+Jp _1 [f-0]“-3d0, i<j  (29) 


Var(u(i)]  =  a-'i)!(n-.i)!-n  /tall2e(i)[e(i)*f]i'1[f-0(i)ln'jde(i) 
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Solution  of  the  Expected  Value  and  Covan  ance  Equations 

Equations  (28) ,  (29) ,  and  (30)  cannot  be  evaluated  in 
terms  of  elementary  functions,  but  they  can  be  evaluated  by 
numerical  integration  techniques  with  the  aid  of  an  electronic 
computer.  The  integral  factor  of  Equation  (28)  does  not  con¬ 
verge  for  values  of  i=l  or  n,  and  Equations  (29)  and  (30) 
do  not  converge  for  values  of  i=l,  2,  n,  and  n-1.  Barnett 
(Ref  1:1209)  points  out  that  for  these  values  of  i,  the 
cos  0  and  cos  0  factors  in  the  denominator  of  the  respective 
integrands  are  dominant  as  0  and  0  tend  toward  their  limit¬ 
ing  values,  and  it  is  therefore  apparent  that  the  means  of 
the  first  and  last  order  statistics  do  not  exist,  and  the 
variances  and  covariances  of  the  first  two  and  last  two 
order  statistics  do  not  exist. 

The  task  of  evaluating  the  remaining  quantities  is  some¬ 
what  reduced  due  to  the  symmetry  of  the  Cauchy  distribution. 
For  the  expected  values 

E[U(i)j  «  -  Efu^i-i))  for  even  n,  ij*l  or  n  (31) 
Elu(n+l)/23  =  0  for  odd  n  (32) 

and  for  the  covariances  of  the  standardized  order  statistics: 

Covlu(i , j) 3  =  Cov[u(n*2-i,n+l-j)]  *  C°v fu( j ,i) 3 

=  Covlu(n+2- j ,n+2-i) 3 

for  i<j  and  i  0  1,2 ,n,  and  n-1  (33) 
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Even  with  this  reduction  in  the  number  of  quantities  to  be 
computed,  the  integrals  must  be  approximated  by  an  elec¬ 
tronic  computer.  Barnett  (Ref  1)  calculated  these  expected 
values  for  n=3(l)20  and  the  variances  and  covariances  for 
n=5(l)20  to  four-decimal-place  accuracy.  Although  the  ex¬ 
pected  values  and  variances  can  be  calculated  quite  effi¬ 
ciently,  the  double  integration  in  the  covariance  expression 
requires  extensive  computer  time.  Stark  (Ref  16:44-48)  cal¬ 
culated  and  tabled  these  quantities  for  sample  sizes  of 
5(1)20  to  six-decimal-place  accuracy.  The  expected  values 
and  covariances  calculated  by  Stark  were  read  into  the  main 
programs  in  Appendix  C  for  the  calculation  of  the  linear 
coefficients  for  conditional  best  linear  invariant  estimation 
of  the  parameters. 
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IV.  I.incar  Parameter  Estimation 


Introduction 

The  purpose  of  this  thesis  is  to  develop  the  coefficients 
required  to  obtain  best  linear  invariant  estimates  of  the 
parameters  of  the  Cauchy  distribution.  This  chapter  explains 
linear  estimation,  develops  the  estimator  with  the  minimum 
mean  square  error  for  both  parameters,  and  describes  the 
method  used  to  solve  the  resulting  simultaneous  equations. 

The  coefficients  are  calculated  for  the  minimally  censored 
sample,  which  consists  of  the  ordered  sample  less  the  four 

i 

extreme  order  statistics  (the  first  two  and  the  last  two), 
and  for  samples  which  are  additionally  censored  from  above 
and  symmetrically.  In  addition  the  values  of  the  mean  square 
error  function  are  computed.  All  of  these  quantities  are 
tables  in  the  Appendices. 

Linear  Estimation 

Linear  estimation  is  a  form  of  estimation  where  the 
ordered  data  are  assigned  weights  and  the  new  values  are 
summed  to  give  the  estimate  of  the  desired  parameter.  These 
weights  are  the  coefficients  and  they  are  calculated  so  as 
to  obtain  the  best  linear  estimator  of  the  parameter.  In 
this  case  the  best  estimator  is  the  one  with  the  minimum 
mean  square  error.  The  expressions  required  to  obtain  these 
best  estimators  are  developed  in  the  next  section. 
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Mean  Square  Error 

Location  Parameter.  From  Equation  (6)  the  conditional 
estimate  of  the  location  parameter  is  given  by 

N 

t* IS  =  Z  aiXi-AS  (34) 

i  =  l 

where  t*  is  the  estimate  of  the  location  parameter 
ai  is  the  coefficient  to  be  computed 
S  is  the  scale  parameter 
A  is  a  constant  to  be  computed 
N  is  the  sample  size 

and  from  Equation  (4)  the  mean  square  error  of  the  estimator 
of  the  location  parameter  is  given  by 

MSEL  =  E[(t*|S-t)2)  (35) 

where  t  is  the  true  value  of  the  location  parameter.  MSEL 
is  the  quantity  to  be  minimized. 

I 

From  Equation  (35) 

MSEL  «=  E[(t*|S)2-2t(t*|S)+t2]  (56) 

From  Ref  10:134 

Var  [xj  *  E [x2 J - (E  [x] ) 2  (37) 

therefore  Var[t*|SJ  =  E((t*IS)2]-(E(t*|S))2  (38) 

and  MSEL  =  Var[t*jSJ  +  (E[t*|S]-t)2  (39) 

Now  by  the  substitution  of  liquation  (34), 
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WSEL  =  V'arM^iX^j-ASl  +  CEt^aiX^j-ASJ-t)2  (40) 

Freund  (Ref  6:173-174)  shows  that  for  a  given  set  of  random 

variables  xx,  x2,...xn  where  Y  =  £  a  linear  combina¬ 

tion  of  N  random  variables. 


E[Y1  “  1f1aiE[xi]  C«) 

and  VarfY]  .  j^Var fx(1) I**t<raI.jCBrIx(i,x{J)  1 

for  i  <j  (42) 

Now  using  relations  (41)  and  (42)  Equation  (40)  becomes 
N 

MSEL  =  if1aiVarU(i)3+2E^aiajCov[x(i)x(j)]  + 

Cif1aiECxCi)3’AS-t)2  C43) 

Hie  expected  values,  variances  and  covariances  developed 
m  Chapter  III  were  for  standardized  order  statistics,  where 

the  standardized  order  statistic  was  defined  by  Equation  (19). 
Therefore 

x(i)  *  SU(i)+t  (44) 

Efx(i)I  =  SElu(i)3+t  (45) 

Varlx(i)]  =  s2Var[u(i))  (46) 

Covtx(i)x(j)3  s  s2Cov(u(i)u(:j))  (47) 

Now  define  the  following  symbols: 

Let  H  *  E[«(i)l 
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a-i  =  Varfu^] 

M  Cov[u(i)u(j)] 

Now  by  substituting  these  relations  into  Equation  (43) 

N  N  N  N 

MSEL  =  Z  a?s2a-,+2  Z  Z  aia.s2an  +  [  Z  a$  (sua  +t) -AS-t] 2  (48) 
i=l  1  11  i=l  j=i+lJ  J  i=l 

N 

and  adding  the  constraint  that  Z  aj=l  and  applying  this 

i=l 

equation  to  a  minimally  censored  Cauchy  sample  of  size  N-4 

(the  two  extreme  order  statistics  are  removed  from  each  end), 

the  following  equation  is  obtained 

N-2  ?  N-2  N-2  N-2 

MSEL  *=  s2(  Z  a|a^^+2  I  Z  aiaj°ij+[  £  aivi"A]^)  (49) 
i=3  i=3  j  =i  +  l  i=3 

Scale  Parameter.  An  expression  similar  to  Equation  (49) 
can  be  developed  for  the  scale  parameter  where  the  mean 
square  error  of  the  estimated  scale  parameter  is  given  by 

MSES  =  E[(S*|t-S)2]  |  (50) 


where  S  is  the  true  value  of  the  scale  parameter, 

S*  is  the  estimate  of  S,  and 
t  is  the  true  value  of  the  location  parameter. 

N 

Let  S*|t  =  Z  djX(i)-Dt  (51) 


and  substitute  into  Equation  (50)  to  obtain 

N  N  , 

MSES  «  Var[  J:^dix(i)-Dt]  +  (EI  ^diX(i)-Dt]-S)2 


(52) 


22 


GAM/MATH/ 72 -3 


And  by  the  use  of  Equations  (42),  (45),  (46),  and  (47), 

Equation  (52)  becomes 

N  N  N  N 

MSES  =  E  d?s2cr.-,.  +  2  E  Es2did-oii+[  E  di (sai+t) -Dt-S] 2  (SL) 

i=l  1  11  i*l  j=i+l  J  J  i=l 

N 

Now  by  adding  the  constraint  that  E  dj  =  D  and  considering 

i=l 

a  minimally  censored  Cauchy  sample  the  mean  square  error  of 

the  estimate  of  the  scale  parameter  is  given  by 

,  N-2  -  N-2  N-2  N-2 

MSES  «  s 2 (  E  dfoii+2  E  E  did^ii+I  E  diPi-1]2)  (54) 

i-3  i=3  j  =i+l  J  J  i=3 

Equations  (54)  and  (49)  are  the  required  expressions  to 
compute  the  mean  square  error,  but  in  this  case,  it  is  de¬ 
sired  to  minimize  these  quantities. 

Minimization  of  the  Mean  Square  Error 

Both  expressions  for  the  mean  square  error  contain  the 
scale  parameter  squared.  At  this  point,  the  problem  is  to 
determine  the  values  of  the  a^’s.  A,  dj's  and  D  which  min¬ 
imize  the  respective  MSE  function,  and  these  values  will  be 
the  same  if  the  functions  are  minimized  without  the  (s2) 
term. 

Taylor  (Ref  17:198)  describes  Lagrange’s  method  of  min¬ 
imizing  a  function  of  several  variables  subject  to  a  con¬ 
straint.  In  applying  this  method  the  original  function  is 
modified  by  adding  the  constraint  equation  multiplied  by  a 
Lagrangian  multiplier,  and  then  taking  partial  derivatives 
of  the  function  with  respect  to  each  variable  and  multiplier. 


23 


GAM/MATil/72-3 


The  resulting  derivatives  are  set  equal  to  zero  to  form  a 
set  of  simultaneous  equations.  These  equations  are  then 
solved  for  the  values  of  the  variables  and  multipliers  which 
minimize  the  function. 

Matrix  Equation.  To  develop  the  matrix  equations  for 

the  calculation  of  the  required  coefficients,  Equation  (49) 

N-2 

is  modified  by  adding  the  constraint,  E  ai=l,  to  give 

i=3 

N-2  ?  N-2  N-2  N-2  , 

L  «*  E  axOii+2  E  E  a^ajaij+I  E  aipj_-A]z+X(Eai-l)  (55) 
i*=3  i=3  j=i+l  J  1=3 

Now  if  N=7  is  the  total  sample  size  and  m=3  is  the  size  of 
the  sample  after  censoring,  application  of  the  Lagrangian 
method  results  in  the  following  set  of  equations. 

3L 

7 X  =  a3+a4+as-l  =  0 
3L 

3ft  a-a3Ji3-a4,«4-a5P5+A  =  0 

2  A 

3ft3e  a3(o33+P3)+a4(c54  +  H3P4)  +  a5(o35  +  {i3P5)-Ap3+-  =  0  (56) 

3L  2  X 

3ft4=  a3(a43+p4y3)  +  a4(a44  +  y4)+a5(o45+p4y5) -AM4+2  85  0 

3  L  2  A 

3ft5c  a3^°S3+^5y3)  +a4  (<J54+»J5i14)+a5  (255+115) -  Ay g+y  =  0 

The  above  equations  in  matrix  form  are: 
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0  0 
0  1 


m»3 


'H 


■**S 


1  -V3  (<y33+*‘3)  (<J34+JJ3JJ4)  (°35+,13M5)  a3  =  0  (57) 


1  -P4  Cct4  3+  ^4 11 3^  C‘r44+lJ4)  (045+^4^5) 

2 

1  -jig  (^53+u5p3)  (054+PSV4)  (o55+p5) 

A  similar  procedure  for  Equation  (54)  results  in  the  follow¬ 
ing  matrix  equation  for  the  coefficients  of  the  scale  par¬ 
ameter. 


tv)  I  >> 

1*  * 

1 

A 

0 

a3 

= 

0 

a4 

0 

a5 
•  « 

0 

.  * 

033+P3)  (034+P3H4) 

2 

(043+M4U3)  (044+114) 

(053+1.5113)  (°S4+ySy4) 

t 

where  D  =  d3  +  d4  +  d3 


C035+P3U5) 

(a45+M4H5) 

2 

(055+P5) 


d3 

< 

^3 

d4 

s 

y4 

« 

d5 

« 

y5. 

(58) 


The  matrix  Equations  (57)  and  (58)  include  matrices  of 
the  expected  values  and  covariance  of  standardized  order 
statistics  and  the  column  vectors  of  the  desired  variables. 
These  equations  can  be  solved  for  the  column  vector  of  m  +  2 
variables  for  the  location  parameter  and  m  variables  for 
the  scale  parameter. 

Solution  of  the  Matrix  Equations.  The  above  matrix 
equations  were  solved  for  sample  sizes  of  N  **  5(1)20  on  the 
CDC  6600  computer.  Basically,  the  main  Fortran  program  reads 
in  the  values  of  the  means,  variances,  and  covariances,  and 
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the  matrix  equations  arc  calculated  and  then  solved  by  a  sub¬ 
routine  to  the  main  Fortran  program.  This  subroutine  is  a 
modification  of  the  "Matrix  Equation  Solver"  Fortran  extended 
subroutine  due  to  the  Computer  Science  Center,  Wright-Patter- 
son  AFB,  Ohio.  The  values  of  the  linear  coefficients  and 
constants  are  tabled  for  each  N  and  M.  In  addition,  the 
value  of  the  mean  square  error  function  is  tabled  for  each 
sample  of  size  N  and  M.  The  values  of  the  mean  square  error 
functions  are  calculated  from  Equations  (49)  and  (54). 

Equation  (49)  can  be  rearranged  to  give: 

N- 2  N-2  N-2 

MSE  «=  E  E  a^ajcr^j+(  E  a^y^-A]^  (59) 

i=3  j  ~3  1=3 

and  a  similar  expression  can  be  developed  for  the  MSE  of  the 
scale  parameter.  With  the  known  values  of  a^  and  A  from  the 
solution  of  Equation  (57) ,  the  MSE  is  easily  computed  in  a 
subroutine  to  the  main  program. 

Additional  Censoring 

Censoring  from  Above.  The  minimally  censored  sample  was 
defined  as  the  basic  ordered  sample  of  size  N  with  the  two 
extreme  order  statistics  censored  from  each  end  of  the  sam¬ 
ple.  Additional  censoring  was  accomplished  by  censoring  the 
sample  from  above  so  that  M  (the  number  of  sample  values 
remaining  after  censoring)  decreases  from  (N-4)  to  (1)  for 
each  sample  size  (N) .  The  value  of  the  mean  square  error 
function  was  also  tabled  for  each  M.  The  values  of  the  co¬ 
efficients  and  MSE  for  the  minimally  censored  sample  and 
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additional  censoring  from  above  are  given  in  Table  I,  Ap¬ 
pendix  A. 

Symmetric  Censoring.  Each  sample  of  size  N=5(l)20  was 
also  censored  symmetrically  from  both  ends.  In  this  case 
M  ranges  from  N-4  to  2  for  even  sample  sizes  and  N-4  to  3 
for  odd  samples.  The  symmetric  censoring  was  terminated  at 
M=3  since  with  M=1  the  estimate  of  the  location  parameter 
is  the  median  and  there  is  no  information  available  to  com¬ 
pute  the  estimate  of  the  scale  parameter. 
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V.  Use  of  the  Tables 

Introduction 

The  results  of  this  thesis  are  Tables  I  and  II  included 
in  the  Appendices.  With  these  tables  the  user  is  able  to 
calculate  the  best  linear  invariant  estimates  of  the  loca¬ 
tion  and  scale  parameters  of  the  Cauchy  distribution.  This 
chapter  explains  these  tables,  gives  the  required  procedure 
to  obtain  the  conditional  estimates  or  simultaneous  estimates, 
and  provides  examples  of  these  calculations. 

Explanation  of  Tables  .1  and  II 

Both  tables  include  the  values  of  the  mean  square  error 
function,  and  the  values  of  the  coefficients  and  constants 
which  are  applied  to  the  sample  data  to  obtain  the  estimates 
of  the  parameters.  These  values  are  tabled  for  each  N 
(sample  size)  and  M  (sample  size  after  censoring).  The  co¬ 
efficients  of  the  order  statistics  required  for  estimation 
of  the  location  parameter  are  listed  in  the  columns  under 
**LOCATION**,  and  the  same  values  for  the  scale  parameter 
are  listed  in  the  same  manner  under  **SCAI.E**.  Table  I  is 
used  for  a  minimally  censored  sample  or  for  a  sample  with 
additional  censoring  from  above,  and  Table  II  is  used  for  a 
sample  which  has  been  censored  symmetrically. 

Estimation  Procedure 

Best  linear  invariant  estimation  of  a  parameter  con¬ 
sists  of  the  following  steps: 
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1.  Obtain  the  sample  data. 

2.  Order  the  data  and  determine  N. 


3.  Minimally  censor  the  two  extreme  sample  values 
from  each  end,  and  additionally  censor  as  desired. 

4.  Determine  N  and  M  and  enter  the  appropriate  table 
to  obtain  the  coefficients  and  constants. 

5.  Multiply  the  sample  values  by  their  respective 
coefficients  and  sum  the  result  to  obtain  the 
estimate  of  the  parameter. 

6.  If  the  sample  was  additionally  censored  from  above 
(Table  I)  conditional  estimation  is  required,  and 
the  appropriate  constant  times  the  Renown  parameter 
must  be  summed  with  the  terms  in  step  5. 


Examples 

No  Additional  Censoring.  As  an  example  of  the  use  of 
the  tables  to  determine  an  estimate  of  the  parameters  of  the 
Cauchy  distribution,  assume  that  the  following  data  are 
known  to  come  from  a  Cauchy  distributed  parent  population. 
The  true  value  of  the  location  parameter  is  8.0  and  the  true 
value  of  the  scale  parameter  is  .5. 


Xj  =  5.689853 

x6 

=  9.222433 

X£  =  7.835235 

x7 

=  8.316519 

X3  =  9.641365 

x8 

=  7.902609 

x4  =  7.201119 

x9 

=18.926210 

x5  =  8.739464 

When  the  data  are  ordered  the  following  order  statistics  are 
obtained: 


29 


GAM/MATH/72-3 


(1) 

tz 

5.689853 

x(6)  =  8.739464 

(2) 

c 

7.201119 

X(7)  -  9.222433 

(3) 

a 

7.835235 

x(8)  =  9;641365 

(4) 

= 

7.902609 

x(9)  =  18.926210 

(5) 

= 

8.316519 

After  these  statistics  are  censored  the  subset  x^  through 
X(7)  remains  with  N=9  and  M=5.  From  Table  I,  the  following 
coefficients  are  obtained: 

a3  =  -.067277  a6  =  .245395 

a4  =  .245395  a7  *  -.067277 

a5  =  .643765  A  =  0 

Now  the  estimate  of  the  location  parameter  t*  is 

=  a3x(3)+a4x(4)+a5x(5)+a6x(6)+a7x(7)  (60) 

And  when  this  calculation  is  carried  out  t*  =  8.290177. 

To  estimate  the  scale  parameter  the  required  coeffic¬ 
ients  aie  obtained  from  the  same  table. 

d3  =  -.153945  d6  =  .369546 

d4  =  -.369546  d7  =  .153945 

d5  «  0  D  =  0 

Now 

s*  »  d3X(3)+d4X(4)+d5x(5)+d6x(6)+d7x(7)  (61) 

and  when  these  calculations  arc  carried  out  s*  =  .522809. 
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It  can  be  seen  by  inspecting  the  tables  that  the  esti¬ 
mation  is  conditional  estimation  when  the  sample  is  censored 
additionally  from  above.  In  all  other  cases  the  values  of 
A  and  D  are  equal  to  zero. 

Additional  Censoring  from  Above.  To  demonstrate  the 
conditional  estimation  procedure,  consider  the  data  from  the 
previous  example.  Additional  censoring  from  above  will  be 
performed  so  that  K- 9  and  M=3.  Now  from  Table  I: 


a3  = 

-.070500 

a5  =  .82517 

a4  = 

.245330 

A  =  -.032654 

d3  = 

-.231235 

ds  =  .779717 

d4  = 

-.588770 

D  =  -.040287 

t*|s 

"  a3x(3)+a4x(4) 

+a5x(5) "As 

(62) 

s*jt 

c  d3x(3)+<*4x(4) 

*'d5xCS)  _1)t 

(63) 

When  these  calculations  are  carried  out,  the  following  esti¬ 
mates  are  obtained:  t*|s  *  8,265232  and  s*|t  =  .342227. 

The  same  procedures  apply  to  the  estimation  of  the  parameters 
with  additional  symmetric  censoring  with  the  coefficients 
from  Table  II.  The  reader  will  notice  that  the  value  of  the 
MSB  of  the  estimate  increases  as  less  information  is  con¬ 
sidered.  For  the  location  parameter  estimator  of  the  last 
example  the  MSB  increases  from  .38655  to  .39752  as  M  de¬ 
creases  from  5  to  3, 
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Simultaneous  Estimation.  Due  to  a  technique  suggested 
by  Herman  for  unbiased  estimation  (Ref  8),  it  is  possible 
to  use  Table  I  for  simultaneous  estimation  of  the  parameters 
when  neither  is  known. 

If  t  «  the  simultaneous  estimate  of  the  location  parameter 
s  **  the  simultaneous  estimate  of  the  scale  parameter 
t*  and  s*  are  as  defined  earlier 


then 


and 


M+2  M+2 

__ 


M+2  M+2 

s  -  ±ls  dix(i)‘Dif3aix(i) 


(64) 

(65) 


As  an  example  of  this  technique,  consider  the  data  from 
the  example  with  "additional  censoring  from  above". 


a3  = 

-.070500 

d3  = 

-.231235 

X(3)  =  7.855235 

a4  c 

.245330 

CU 

-b. 

ii 

-.588770 

x(4)  =  7'.  902609 

a5  a 

.82517 

*5  * 

.779717 

x^5j  =  8.316519 

A  = 

-.032654 

D  = 

-.040287 

% 


-  8.248905  +  .000651 

r  e  .998684 


=  8.200428 


r  _  .019931  +  .332324 

T9DW84 


.352719 


In  summary,  it  can  be  seen  that  although  this  technique 
of  estimation  is  conditional  estimation,  the  conditional  re¬ 
quirement  is  only  necessary  when  estimating  from  sample  data 
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that  have  been  censored  additionally  from  above.  The  method 
provides  a  simple  and  efficient  procedure  to  estimate  the 
Cauchy  parameters. 


33 


GAM/MATH/ 72 -3 


VI .  Summary 

The  objective  of  this  thesis  was  to  develop  a  table  of 
linear  coefficients  which  could  be  easily  applied  to  sample 
data  to  obtain  the  conditional  best  linear  invariant  esti¬ 
mates  of  the  location  and  scale  parameters  of  the  Cauchy 
distribution.  These  estimates  are  best  in  the  sense  that 
they  process  minimum  mean  square  error.  The  coefficients 
and  constants  required  to  calculate  these  estimates  for 
minimally  censored  samples  and  samples  with  additional 
censoring  from  above  are  tabled  in  Appendix  A.  The  same 
values  for  samples  with  additional  symmetric  censor’ ng  from 
above  and  below  are  tabled  in  Appendix  B.  In  addition  the 
.  value  of  the  mean  square  error  function  for  each  sample  size 
is  also  included  in  these  tables.  The  computer  programs 
required  to  calculate  and  table  the  above  data  are  included 
in  Appendix  C. 

This  paper  also  includes  a  brief  review  of  the  previous 
work  in  estimation  of  the  parameters  of  the  Cauchy  distribu¬ 
tion  from  1961  through  1970.  Much  of  this  effort  was  ac¬ 
complished  at  the  Air  Force  Institute,  of  Technology  under 
the  direction  of  Professor  Albert  H.  Moore  and  sponsorship 
of  Dr.  H.  Leon  Harter.  These  works  are  concerned  primarily 
with  unbiased  estimates  of  the  parameters. 

The  Cauchy  distribution  and  its  peculiarities  are  dis¬ 
cussed  with  a  review  of  several  ways  in  which  the  distribu¬ 
tion  may  be  generated.  The  order  statistic  theory  required 


34 


GAM/MATH/ 72 -3 


to  obtain  the  expected  values,  variance,  and  covariances 
of  the  order  statistics  is  also  reviewed.  The  mean  square 
error  function  for  the  estimate  of  the  parameters  is  de¬ 
veloped  and  minimized  by  Lagrangian  techniques  to  obtain  the 
matrix  equations  required  to  calculate  the  linear  coeffic¬ 
ients. 

The  report  is  concluded  with  an  explanation  of  the 
tables  and  several  examples  of  the  application  of  these  co¬ 
efficients.  A  technique  of  simultaneous  estimation  of  both 
parameters  of  the  distribution  is  also  presented  with  an 
example  of  the  technique. 

The  method  of  estimation  presented  in  this  paper  and 
the  attached  tables  of  linear  coefficients  provide  a  simple 
and  efficient  method  of  obtaining  either  conditional  or 
simultaneous  best  linear  invariant  estimates  of  the  param¬ 
eters  of  the  Cauchy  distribution. 


\ 
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TA^t r  I 

f'OcFF ICIr. NTS  FOR  3£ST  CCNTITTCNAL  F*TI«ATICN  CF  THF 
LOCATION  AND  ?CALr  FADAweTSFc'  OF  TH£  CAUCHY  CISTRI^UTICN 
(WITH  ADOll  ICNAL  CFNSC-'IMR  CC,CH  AprvF) 
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TA^L:  I 

COrFriciF‘ITS  FC7  FFST  CCNdTICNH  rSTTMCTICK  OF 
LOCATION  A NC  SCALE  F  A  ^ Aw  H  7  rRS  OF  THF  CALCHY  CIST  9T*hjti  r*.; 
(WITH  AUCTTICKAL  C2NFC3ING  F^Ct'  APCVn 


**  LOCATION  **  **  SCALE  ** 
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Tft^LF  I 

rcrFFTCxrrnc  fc?  erFt  ^csfittcnal  fstjmattcn  cf  th^* 

L CCA T ION  ANH  <C AL *=•  CARAvET^rS  CF  TH ^  CflCMY  CIST9InUTICN 
(WITH  APOIlICNfL  CF NrC9  TNG  FRC»'  ArCVf.) 
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T-^LF  I 
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Pf  ICTr NTC 
TON  A^H  «: 
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,1r ST  CCKPTTTCUH 
F  A  R  A  X  £  T F  S  CF  Tn 
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ESTIKA-ICN 

*■  calchy  r.ir 

FkOH  APCVE) 
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T  r*Ir!L‘TT  CN 
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fa°LE 

FOR  t?FcT  CCf.TITKNAL  FSTIRATICN  CF  th<* 
LOCATION  A,:D  T^ALt  cA'RAy£TFPS  OF  TFT  CALOFV  CITTRIPUTICN 
{ VI TIM  AOQITI.OKAL  OENFCP ING  fPOK  A<?CVn 


**  LOCATICN  **  **  fCAL-  ** 

N  M  M?F  I  rCSP.  M?r  I 
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TArLF  I 

COr-FcI'!i;  NTS  FO?  Pf'JT  OCNriTTf'KAL  L^TiyflTICN  r*  T~- 
LC9AT1CN  At.:?  £  C  A  L  c  PARAMSTee?  OF  ThF  CALOwy  ri<Ts:“'J*IC*. 
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T  A  nL  E  I 


HOrFF  IDIr  NTC  FQ°  PE<'T  CCNTITICMAL  ESTIMATION  CF  THE 
LC^ATICN  AND  SCALC  FA9AM°Tt.F5  CF  THE  CAL’CFY  CISTRIStJTICN 
(WITH  ADDITIONAL  C-NSOFING  FrOH  A0CVF) 


**  LOCATION  **  **  SCALE  ** 


H 

M$E 

I 

COSF. 

MFE 

I 

°CEF. 

15 

3 

.4764? 

■? 

-.  36?:-6F 

.35032 

•J 

V- 

- • P45  F*  9 

4 

-.13346? 

4 

-  *  1 894"J5 

5 

1. 2*5  25t7 

5 

-.  315  ?: 3 

A 

-.5722°? 

C 

-.553856 

15 

? 

l.i?32K 

3 

-« 174  992 

.41034 

3 

-.046684 

4 

1.1749°? 

4 

-.3?24.  2 

A 

-1.1574'-? 

C 

-•419rfi  6 

15 

1 

5.24717 

3 

1 .5  j;  ,cr 

• 5244u 

V 

-.218*17 

A 

-2. 1*1462 

c 

-.218:3? 

IF 

12 

.i?7i9 

*7 

-.3? *14? 

•15819 

7 

-.014524 

4 

-.337.67 

4 

„•  6  0  4  7  7 

c 

+ 

- , C 1 35° 1 

5 

-.1313*7 

e 

,3??c6° 

6 

-.186477 

7 

,2  *0  927 

7 

-.17221° 

* 

.2  =  Fcit. 

a 

-. 07? 7° 8 

9 

,29651 4 

c 

,->•»?  7c« 

1C 

.2:0627 

13 

.172519 

11 

, 173^6  ? 

11 

.1*647? 

12 

-. J13597 

12 

•1’IJF 7 

1  0 

-.317.67 

1? 

.«6?47? 

14 

32514? 

14 

.  i>14c?4 

A 

•  o«3  :tc 

c 

-•oocc-.c 

16 

11 

•  1 6  .*  3  7 

*7 

-«12C5°f 

.  15885 

-.  0146*6 

4 

-.337344 

4 

-.06 J722 

5 

-.013674 

5 

-.131=35 

6 

•374?*? 

6 

-.  1.373:  i 

7 

.212274 

7 

-,173rl° 

•3 

.2  9*96.4 

8 

-•  w71 22? 

9 

.29” :?? 

e 

* 

.  c7???: 

1  ) 

.232474 

13 

.172721 

1 1 

•374cei 

11 

.1*7  O  '  5 

12 

-  •  .0  1  7  7  2  *• 

32 

.  l?le:.0 

13 

-•06FC71 

1  ’ 

.012*11 

A 

.  3  *,5*4? 

C 

.*35*12 

■  ',1 


1 


GW/MftTM/72-3 


T  A  nLF  T 

CCfrFriCIrHT?  FG9  95ST  CCNdTICNAL  FSTIMM  ICN  CF  THE 
LOCATION  AND  °CAL5  FA9AM2TFPS  OF  THE  CAUCHY  CTST9IPUTICN 

(with  additional  censoring  f°oh  aocvf) 
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ta  fl  r  I 

nC?FPICIrNTC  FOR  "f  ST  CCNTJT  TONAL  FSTI^ATICN  CF  Tf-F 
LCCAT ION  AND  SCAl"  FARAMcTpPS  ?F  THC  CAICHY  CIST°I"LTTCN 
(WITH  AOOITIONAL  CLOSING  cP0K  A9C VE ) 


**  LOCATION  ** 

N  M  MFF  I  "OF" 


**  SCALE  ** 

MFC  I  COFF. 


16 

7 

.17397 

3 

-.32090? 

• 235"1 

3 

-. 02195? 

A 

-.030640 

A 

-. J91766 

3 

-.014272 

5 

-.2cC7C5 

A 

«  U76  A7‘! 

6 

-.20‘*77-’ 

7 

.20 36"? 

7 

-.275791 

3 

.2C842? 

8 

-,132e?9 

9 

.430397 

9 

1.0A1M6 

A 

•  32 1 '1 6 

C 

.020499 

16 

6 

.1«459 

3 

-.022292 

• 2°46 " 

7 

V 

-.027751 

A 

-.341601 

A 

-.116240 

5 

-.017  IFF 

5 

-.255700 

6 

6 

-•  271491 

7 

,214  j£o 

7 

-.  ?62p46 

3 

,790269 

8 

1.014<R6 

A 

-.  27  5  23? 

C 

-,12'C40 

16 

5 

.21736 

3 

-."2675l 

.74088 

3 

-.  37?,56 

A 

-.  351C.?f 

A 

-•  14"  c66 

5 

-. 025724 

5 

-.31jr.?6 

6 

.079:02 

6 

-. 454714 

7 

1.324427 

7 

.567076 

A 

-.221221 

C 

-.371205 

16 

A 

.  3  r  0  4  c. 

3 

-. j  37576 

.36833 

7 

V 

-.035637 

A 

-.  J75517 

A 

-.15C2<2 

6 

— .04925" 

5 

-.333147 

6 

1.1620F2 

6 

-.C’403? 

A 

-.45163C 

P 

-•  5536?9 

16 

3 

.92126 

7 

-  ,3  r.  q  f.  t.  4 

. 36844 

3 

-.035661 

A 

-.I4475r 

A 

-, 153 "28 

5 

1.213' 28 

5 

-, 25555  C 

A 

-.765624 

C 

-.541179 

1* 

2 

1.27366 

3 

-•177.. 24 

.3969a 

7 

-.037151 

A 

1.177  .24 

A 

75°?54 

A 

-1 . 2  7  2  lr 8 

P 

-.7965: 4 

GflfVMATW/72-3 


TA*L£ 


COEFFICIENTS  FC9  CEST  CCNTITIOf'H  ESTIMATION  CF  THE 
LOCATION  AMD  SCALE  FA<?Am£TFP5  OF  THE  CAUCHY  CIST9IrUTICN 
(WITH  AOOITICNAL  CENSORING  F^CM  ABC VE ) 


N  M 


16  1 

17  13 


17  12 


**  LOCATION  **  **  SCALE  ** 


MSF 

I 

COcP. 

MSE 

I 

CCEF  • 

5.98125 

3 

1.30331:. 

.51914 

3 

-.274294 

A 

-2.353784 

D 

-.234204 

.15440 

3 

-.316948 

.14677 

? 

v 

-.011248 

4 

-.032775 

4 

-•04787® 

5 

- • 372  3F 1 

5 

-,ir7651 

6 

.01927? 

6 

-.  162*19 

7 

.143717 

7 

-• 175 lr  1 

3 

.245754 

3 

-.115015 

q 

.2-3*461 

o 

.  00  0  Cl  3 

1G 

.245754 

10 

.  115  ri5 

11 

.141717 

11 

.175151 

12 

, 3  ?c  272 

12 

•  1^3  51  9 

13 

- . ?  2c "61 

12 

.137*51 

14 

- .  3 1 7  77r 

14 

.047578 

i5 

-.315945 

3.5 

.011248 

A 

♦  ::kgc 

C 

. *32  C9  rt 

.15543 

3 

-.317rfC 

.14723 

7 

V 

-.311254 

4 

-•3** **995 

4 

-. *43930 

5 

-.3224°r- 

5 

-.137505 

6 

•  3  29 66  4 

6 

-.16445® 

7 

.144729 

7 

-.175736 

3 

.  ?  4  7  4  f  1 

8 

-.115439 

9 

,293577 

0 

a# 

-.  CCG  C°6 

13 

,247572 

10 

.115247 

11 

.144  ice 

11 

.176543 

12 

, 329592 

12 

.164257 

12 

-.02273° 

12 

.107606 

14 

-.0  5  nee'7 

14 

•  064698 

A 

. J'46t* 

C 

.  CC4422 

GAfVMATH/72-? 


TAPIE  T 

COEFFICIENT ?  PCS  pp$T  CONDITIONAL  ESTIMATION  CP  THE 
L CCA TIC N  ANO  SCALC  PARAMETERS  CP  TFE  CALCHY  CTSTRIPL'TICN 
(WITH  AnniTIONAL  CENSORING  FROM  AnCVE) 


N  M 


17  11 


17  10 


17  9 


**  LCCATICN  ♦*  **  SCALE  ** 


MC£ 

I 

COEP. 

MSE 

I 

CCEF. 

1674? 

*7 

-•017275 

.14903 

3 

-.311572 

4 

-.93441  4 

4 

-  •  3  4  3  P  3  0 

5 

-.322725 

c 

* 

-• ln9Pll 

6 

.  •)41'2fi7 

6 

-.167263 

7 

.146?f2 

7 

-• 178  °46 

d 

•  251 »C  ? 

e 

-.117713 

Q 

.294^74 

9 

- .  0  0  y  F  7  3 

10 

.251412 

1C 

.116287 

11 

•  1 4  7  5  A  ’ 

11 

.177536 

12 

.041261 

12 

.165569 

13 

-.393*8? 

13 

.103291 

A 

.01954/, 

D 

.013577 

ir9?r. 

3 

-.017465 

•  1667? 

■j 

\. 

-.012134 

A 

-.02477? 

4 

-.  C51213 

c 

* 

5 

-.115222 

D 

.0409C9 

6 

-.1756 r 3 

7 

.14? £93 

7 

-.1332C1 

3 

.254536- 

8 

-.124639 

9 

.2960  Pc 

c 

-•  C02F56 

13 

.255471 

13 

.119466 

11 

.1506?? 

11 

.18311 4 

12 

-.074439 

12 

.439684 

A 

.043159 

C 

.042490 

1596.fi 

7 

-.317512 

.17250 

3 

-.013738 

A 

-.32484? 

4 

-.056545 

5 

-.022872 

5 

-.127766 

6 

.341,fC 

6 

-.194567 

7 

.145717 

7 

-• 2J°4?1 

3 

•°559r2 

e 

-.140655 

9 

.3^0799 

9 

-.007614 

10 

.257266 

10 

.1253-2 

11 

.;  73257 

11 

.63968? 

A 

,05594? 

c 

.064757 

r>4 


.71695 


GAN/WATH/72 


tahle;  i 

COEFFICIENTS  50°  6r?T  CCNPIT  T CMAL  ESTIMATION  CF  THE 
LCCATIOM  AND  SCALE  FAFA^TEFS  CF  THE  CAUCHY  C I STRI6UTT CN 
(WITH  ADDITIONAL  CENSORING  FROM  ArCVE ) 


LOCATION  **  **  SCALE  ** 


N 

M 

njcr 

I 

CCEF « 

USE 

I 

CCfF. 

17 

3 

,57-'1«o 

3 

-.068236 

, 35  C  46 

3 

-.026309 

4 

-.15273? 

4 

-.1235*6 

5 

1.2216? 

5 

-.376428 

A 

-.655725 

C 

-.525714 

17 

2 

1.42468 

3 

-.17*705 

.38630 

■J 

V 

-.G30C92 

4 

1.17*715 

4 

-.3453? 1 

A 

-1.3844 Cl 

D 

-.3*539? 

17 

1 

6.7f6*P 

3 

1 . 0  2  C  v  C  0 

.51497 

-.192132 

A 

-2.6244EC 

r 

-.19213? 

13 

14 

.1434Q 

-.13427? 

.13667 

3 

-.00899? 

4 

-«92G3r9 

4 

-.628353 

5 

-.''256C5 

c 

-.083544 

6 

.  3 1 6 «;  7  3 

6 

-.142  01-: 

7 

.  5  9  c  ■»  3  r 

7 

-.157454 

3 

.1S57?P 

f 

-.133745 

o 

. 26  0  ’74 

c 

* 

-•  '’63  68  4 

10 

.26^334 

10 

.  0538*4 

11 

.19533* 

11 

•  13  0  345 

12 

.?9°335 

12 

.167494 

4  *? 

■a.  V 

,  C 1 6  2  7  2 

13 

•  14?  r 1 0 

14 

-.02661 o 

14 

.033544 

15 

-. 33  a  ;ca 

15 

.  «•  3  3  3  5  7 

16 

-.01427? 

16 

.  0CS997 

A 

1 

P 

-  •  C  •?  -  C  f  « 

GAh/M ATH/72 


7  A  9LF  T 

COEFF  JCTrNT ^  FOR  BEST  CONOITICNAL  ESTTMATICN  Cc  THE 
LCCATTCN  ANO  SCALE  PA  s/WETERS  OF  TFE  CAUCHY  C I 57RInUTI CU 
(WITH  AODTT ION AL  C£NcOnING  F»OM  APCVE) 


N  M 


18  13 


13  12 


**  LOCAT  IO‘J  **  **  SC  ALC  ** 


M?r 

I 

OQ£F. 

KSF 

I 

CCEF. 

1443? 

7 

-.ni4sr4 

.13719 

3 

-.oC°"19 

4 

-•97Cn71 

4 

-.938445 

5 

T267TF 

5 

-.038757 

6 

.916464 

6 

-.  1423t>K 

7 

.OCC04- 

7 

-.167969 

3 

♦  IBF'7!? 

8 

-.138706 

q 

.261 PIS 

c 

-.954*64 

i  j 

.25194* 

13 

.053528 

n 

.196*05 

11 

. 1?  Bc  ft  C 

12 

•  i:-3r7P 

12 

.1677»1 

13 

•  u  1 n  5  59 

12 

.14^224 

14 

-.12* US 

14 

.t6«6C9 

13 

-.u"laP2 

15 

.  j  5 1  p  8  9 

A 

.  0  :  *»  7  7 1 

C 

.o03585 

14602 

3 

-.014f2? 

.13639 

3 

-.609132 

4 

-.033923 

4 

-.o38529 

5 

-.C27942 

5 

-  .  0  8  9  f  8  3 

6 

♦916744 

e 

-.144180 

7 

•191247 

7 

-.170113 

3 

.199:46 

8 

-• 14" f 37 

g 

.265**24 

5 

0*5*52 

10 

.2654  58 

10 

.054184 

11 

.  1994F'* 

11 

.139761- 

12 

.101366 

12 

.169249 

13 

.017924 

13 

.142315 

14 

-.994077 

14 

.156616 

A 

•016197 

C 

.oict:6 

67 


. . >  vw  *  ■■  ■  l  ^ "  m.  ',  'u  iwwm  1 1  ut ;  i" ,' v»-y X’-'W  V"  “ ™ 


GAM/MATH/72-3 


TAPLP  I 

HOFFP ICIFNTC  FCW  GOST  CONDITIONAL  ESTIMATION  CF  THE 
LOG  A  TICK'  A  md  SCALE  FARAVSTEPS  CF  THE  CAUCf-Y  CIST9InUTICN 
(WITH  ADDITIONAL  CENSORING  PROM  APCVE) 


N  M 


13  11 


is  in 


13  9 


**  LOCATION  **  **  SCALE  ** 


MFF 

T 

COEF. 

MSC 

I 

CCEF . 

14779 

7 

-.314797 

.14293 

3 

-.009470 

4 

J 7 1279 

4 

-.04*276 

5 

-.327310 

5 

-.093249 

6 

.017  ->66 

6 

-. 1496^2 

7 

•  1 32  74C 

7 

-.174/40 

3 

.231904 

8 

-.14649? 

9 

.269173 

c 

-.  C5312u 

10 

.269124 

10 

.054714 

11 

•  2'‘?*PB 

11 

.14211 3 

12 

*134254 

12 

.173431 

l* 

-.024180 

1? 

. 33^527 

A 

♦  337675 

r 

.336691 

1496C 

3 

-.314*77 

.15519 

2 

-.  01C??3 

4 

-.031420 

4 

-.04*2*6 

5 

-  •  02 7  ■'87 

5 

-.1':  0623 

6 

.ir,.?46 

6 

-.162C51 

7 

.13201? 

7 

-•1°1764 

3 

.20391? 

8 

-.1*9293 

9 

.27137? 

c 

-• 065  4ft 2 

1C 

.271929 

10 

. J5F  29  ? 

11 

.235102 

11 

.150025 

1? 

.303*15 

12 

.589066 

A 

.05*111 

C 

.  uG-^eag 

i  a  «  g  r 

3 

-.914*73 

.17627 

3 

-.•J116r5 

4 

-•OM  42F 

4 

-.049770 

5 

-.337263 

5 

-.115169 

6 

.017744 

6 

-.135479 

7 

.1 j36Cfi 

7 

-•  2?’’  S'7? 

3 

•  2  0  3  c  1 5 

8 

-.186c?t, 

9 

.271  ;6* 

9 

0  8 " 4  4  6 

n 

.271919 

14 

•  055  266 

n 

.935945 

11 

.862464 

A 

.^57945 

C 

.063736 

6P 


I'Pf  IH'HlIJMj 


G  Ah'/ M  ATM/ 72 -7 


TA^LE  I 

CCHFcIGIt  NTC  PC°  CEST  CCNFITICMAL  FSTIMATICN  CF 
LCCATICN  A r:0  Sr-ALE  ?A3A”-t?fc  of  Thf  CAUCHY  r ISTRICUTI CfJ 
(WITH  atcitional  CENSORING  FRCK  A5CVF) 


**  LOCATION  **  **  SCALE  ** 


N 

M 

NFF 

I 

rOEF  • 

MSF 

I 

CCEF . 

13 

3 

.15:20 

3 

-.019  IF? 

•  21065 

3 

-•  314  ;.2? 

4 

-.031 66  p 

4 

-.059926 

5 

-.o '7=e? 

5 

-.136972 

6 

.015 9*5 

6 

-.224587 

7 

.103633 

7 

-.206915 

3 

•  2 ' 4 f  34 

8 

-.27C559 

c 

.  2  7  2 1C  1 

c 

-.107163 

10 

.47?  37 « 

1  u 

1,071 7r9 

A 

.019166. 

C 

.  020698 

19 

7 

.  1 5  8  3  G 

3 

-♦C15923 

.25751 

1 

V. 

-,L17?°6 

4 

- ...  3  3  6  6  ■» 

4 

-.374^-3 

5 

-...7  0  756 

c 

-.172223 

6 

.  j  ?.***■» 

6 

-.27‘'-f  « 

7 

.109761 

7 

- ,  3  ?  7  7  8  4 

3 

,209 73 5 

8 

-.29556? 

9 

«7487°4 

o 

1.  J733-2 

A 

-.069329 

D 

-, 106  2°  3 

13 

6 

» 1 6  j  j  4 

7 

-.  j 1326 7 

.30755 

3 

-. u23  C16 

4 

- , 3  39  3C  7 

4 

-.08*445 

<7 

-. 3  36  /A ^ 

5 

-.23 6 2° 7 

G 

.9 13 399 

6 

-.330640 

7 

,  1 1 2  ‘  7  r 

7 

-• 41o  33  4 

3 

.96  a  S8Q 

8 

.733574 

A 

-.  19.4  961 

C 

-.328739 

18 

5 

.?2°2 3 

3 

-,^2361 4 

.72746 

3 

-.022468 

4 

-.051 ??5 

4 

-.  096645 

5 

-.451766 

c 

-.229727 

3 

.  )  ;  C  e  F  £ 

6 

36°?67 

7 

1 . 12',r2p 

7 

,  1 p  4  ^99 

A 

-.371132 

r 

-• 57 C  r2° 

13 

4 

1 

w 

-.9359C9 

.3299? 

V 

-.C227C  2 

4 

-  •  ■}  S 1 '  i-9 

4 

- »  0  9  7  p  0  3 

9 

—  ,  a  c  1  2  c 

E 

-,  228-7°2 

6 

i.?:c  .p-r 

6 

-.236768 

A 

-.^77  ;  f-  N 

0 

-.587454 

co 


6AK/MATH/72-? 


7  A  RL  £  I 

nOEFFI^IFNTS  c03  Rf.ST  OCM'ITICNAL  ESTIMATION  Cp  T«E 
LCCAT10N  AMO  STALE  PARAMETERS  OF  TF  r  CAUCFY  CIMrflPt'TICN 
(WITH  A finiT  IONAL  Crf'SOOIMG  FRCF  A FOVF  ) 


•  •  • 

N 

•  •  •  •  l 

M 

♦  ♦ 

MST 

LOCATION'  ** 

I  OC"F» 

FSF 

SCALE  ** 

I  CCRP. 

13 

3 

.62501 

t 

-.'6.5111 

.72670 

3 

-.022557 

4 

-.15325° 

4 

-.09*28* 

5 

1.227370 

K 

.-.355  44  7 

A 

-.943231 

P 

-• 5J6E21 

1*3 

2 

1.5864? 

7 

-.ncicc 

.37773 

7 

-. 024743 

4 

l.lOClf F 

4 

-. 331158 

A 

-1.494767 

r 

-.355936 

13 

1 

7.6C217 

7 

.51163 

7 

V 

-.181294 

A 

-2.69??9° 

C 

-. 181 2°4 

10 

15 

•12431 

7 

N* 

-  •  u  1  ?  2  7  P 

• 12P21 

7 

v 

-.  0  j 7 2 2 7 

4 

-.:27i°4 

4 

-.031 

5 

- .  j  2  ?. !  1- * 

c 

-.t7?2"8 

6 

,  2  3 1  ?  f  C 

e. 

-.122152 

7 

.  .3  55/' 7 

7 

-.154747 

•3 

,(rrrc; 

e 

-.147151 

0 

,227'92 

c 

* 

-.  ;<r4r.5 

1C: 

.252FT7 

10 

.  ccr-Ou: 

11 

,?23e°? 

11 

. 090459 

12 

.15'5°f 

12 

.147151 

17 

,  0  e  5  5  <•  3 

13 

. 154747 

14 

.  o :  1 2 1  a 

14 

.12'  OF  2 

15 

-  •  ..  c  ?  i  r  c 

15 

.  j?  3  7r  3 

16 

-. j2?1° 4 

16 

.  uil  <*59 

1  7 

-.712270 

17 

•  C  C  7  2  ?  7 

A 

.  w  1  o 

C 

.307000 

COEFFI^If.  MTS  FOR  If  ST  CCNXITIONAL  rSTIMATICN  CF  THF 
LCCATICN  Arn  SCALE  PA^Af'ETFRS  CF  TF5  CAUCHY  C I?T 9 IBUTI CN 
(WITH  ACOITIONAL  CENTRING  FRCF  A^CVE) 


**  LOCATION  **  **  cCf LE  ** 


N 

V 

M<;r 

V 

I 

CCFF. 

MSF 

1 

ccef. 

t 

;  i9 

14 

*2 

w 

-.  Jl22?° 

. 12044 

3 

-.  397241 

r 

4 

-.02?  ?29 

4 

- . u31 116 

[ 

9 

-.32*324 

5 

-.  u73247 

6 

.i,11'52? 

6 

-. 1223P9 

7 

•  C66  la? 

7 

-,155f 23 

3 

.151377 

e 

-.147477 

9 

.224 *L2 

c 

-. j?oer? 

11 

.2f29co 

10 

-.  j09f.-42 

11 

.2949(2 

11 

.  '.93  5c  8 

1’ 

•  151457 

12 

.147259 

1’ 

.160  213 

13 

•154953 

14 

.',"i?r7 

14 

.122 ’5 3 

: 

15 

-.92617? 

15 

.  i*  7  ?  ?  5  2 

IS 

-.04540*5 

16 

.041774 

A 

.or*»:77 

C 

.332924 

19 

13 

« 12613 

2 

-.912471 

.12957 

V 

-.0-7311 

4 

-.327518 

4 

-.C31417 

5 

5 

-.C'74'-56 

o 

,  TO  1273 

F 

-.122538 

7 

. 3669P? 

7 

-•  1565S4 

S 

.1.52177 

8 

-.14»?t;g 

9 

.22747? 

9 

-.  091592 

10 

.2359?? 

10 

- »  00  C  291 

11 

»227f 73 

11 

.001112 

12 

.15351-9 

12 

.14*275 

1? 

.0674P8 

12 

.15599? 

14 

.u3i-af  5 

14 

.1229? 3 

15 

-  • n  9  7  7 1  r-> 

If 

♦12*271 

A 

.0124P3 

r 

•  012  P4  7 

<•1 

GAF/NATH/72-7 


T A nLE  I 

CCEFrinirHTS  FOR  FF5T  rCNDITTCNAl  ESTIMATION  TF  THE 
LOCATION  AND  FPfli£  F A9AM5T  5F  S  OF  ThF  OAfCFV  CIFTRinUTICN 
(WITH  AfJTITICNAL  CENSORING  FROM  A«CVE) 


•  •  t 

N 

M 

**  LO 
MFE 

CATION  ** 

I  COEF. 

** 

MSF 

CCAL5  ** 

I  COST. 

19 

12 

•17781 

3 

-.01 2*22 

.12232 

T 

-.037521 

4 

-.C27E47 

4 

-.  C32223 

5 

-  •  C  “*  8  9?° 

5 

-.676106 

6 

.131376 

6 

-.127198 

• 

7 

.v»67577 

7 

-.161236 

9 

.1 3Eirf 

8 

-.153566 

q 

,2307iC 

c 

-.004606 

1J 

.263732 

3  G 

-•«,311?? 

11 

.221101 

33 

.092629 

12 

.156"?° 

32 

.  151 3 1 6 

13 

.35925* 

12 

.16o;.i 

m 

- .  1 "  2  2f  ? 

34 

,  2 c  ?  6  5  ? 

A 

.032612 

C 

•  «  7 1 6  5  r 

19 

11 

.17SP5 

-.9127J./. 

.14151 

s* 

-.j"75Pi 

4 

-.  3  ?«!/••  0 

4 

-.  1*34254 

5 

-.029 . 9 k 

5 

-.001'?? 

6 

.33156? 

6 

- • 135  2F  2 

7 

.0  2373  p 

7 

-  171709 

6 

•  1  5  c  pz  7 

8 

-. 16437; 

q 

•  2  3  2  9  8  4 

9 

-.132200 

1.1 

.262'.76 

1C- 

-.Dd7232 

11 

.232379 

11 

.095  7,", 

12 

.153459 

12 

.167572 

13 

-.045??- 

12 

.502121 

A 

,054100 

0 

. 355395 

19 

in 

.i:-:enr 

3 

-•U 12 723 

.15694 

w 

-. -98885 

4 

-.329154 

4 

-.036217 

5 

-•■325091 

c 

-. 3  °  0  2 : 9 

6 

•  30197'.. 

6 

-.150  957 

,96  81'F  9 

7 

-.191991 

4 

.157172 

n 

-.194234 

9 

•  2  3  3  4f  c 

c 

-.116933 

n 

.  264  7 

1C 

-  .  ;  ?  8  6  r  q 

n 

•  2  24  4p6 

31 

.103*24 

12 

.3 19?r6 

32 

.760415 

A 

•  02  7  If  6 

n 

.  07126  3 

GAN/H1TH/72-? 


TA^LC  I 

OOEFFICTtMTS  P0R  n£ST  CCNFITICMAL  FSTIMATICN  CF  THF 
LOCATION  AND  SCALE  FA='Af'ETFRS  CF  TFF  CALC^Y  ncTRI?UTICN 
(WITM  ADDITIONAL  CENSORING  FR"K  AnCVF ) 


N  « 

1°  9 


19  « 

19  7 

19  A 


**  LOCATION  **  **  SCALE  ** 


«SE 

I 

CCEF. 

MSF 

I 

CCEF  • 

13943 

*? 

-.032771 

.18263 

2 

-•C102°6 

4 

-  •  0  2  ?  979 

4 

-.044712 

5 

J29276 

5 

-.106600 

6 

.091379 

6 

-. 17?  169 

7 

.058766 

7 

-.22623  6 

8 

•  1 5?13  4 

8 

-.219197 

a 

•  22  24?3 

o 

-.142p63 

l'J 

.262«Ni 

3  0 

-•  91 A  3 . 4 

11 

.245812 

11 

l.CO'1462 

A 

.  342918 

C 

. C 56216 

1“‘J2F 

7 

-.ul3l<‘7 

.22032 

n 

-.n’AC? 

4 

-.029182 

4 

-• «64  231 

9 

-.03  "<.6? 

5 

- .  3  2  ;5  6  7  1 

4 

.00r»£6 

6 

-.216V9 

7 

.C69114 

7 

-. 277 cr ,  * 

8 

.lAO-.tc: 

8 

-.2?224? 

9 

.236498 

c 

-.132250 

1? 

.637C79 

1C 

1. 120  2?7 

A 

-.01641 0 

C 

*-.•325214 

iK5C  C 

3 

- . 31 4 2° 2 

.26445 

3 

-. j15  262 

4 

-.32193  3 

4 

-.065967 

5 

-*  f*23°?i 

5 

-.15658  "* 

6 

-.3C14g6 

6 

-.264T- 

7 

•  373  r19 

7 

-.341 'r7 

8 

•166276 

8 

-•  373993 

9 

.344226 

q 

.934616 

A 

-.116142 

C 

-.1931*3 

,1819? 

3 

-.OK-932 

.29991 

V 

-.0l7p5  1 

H 

-.C28227 

4 

-.  i»75C61 

5 

-.C 41 993 

a 

-.17999? 

6 

-.coe7?^ 

ft 

-.304879 

7 

.3’*575? 

7 

-.396482 

a 

1.329336 

8 

•5537c? 

A 

-.255296 

0 

-.420763 

GAV/^ATW/72-7 


T  A  '’LF  I 


OGl'FF  ICTr  NTS  FCS  CCNFITICNf L  ESTIMAT ICN  Cc  THP 

LCOATIOM  Af:n  SPAL^  FARARFTFF?  CF  THF  CAICMY  C151  R1?UTI CN 
(WJTH  AHniTIOKAL  CZhFOFING  FRCP  ACCVF) 


**  LOCATION  ** 

N  **  MFF  I  fC  =  F 


**  SCALE  ** 

KSE  I  CCEF  • 


19 

5 

.27903 

3 

-.0225*1 

.3121/7 

3 

-.018412 

4 

-• J?1 °64 

4 

-• Q796P4 

5 

-.366761 

5 

-•  139  FF7 

6 

-.017567 

6 

- • 3?22c9 

7 

1.152112 

7 

.026786 

A 

-.437972 

C 

-.573527 

19 

4 

."6471 

■j 

w 

-.025512 

.31216 

7 

V 

-.biSi?? 

4 

-.  J  45*15 

4 

-•  u79783 

9 

-.ni3c« 

c 

-.190121 

6 

1.219:2  4 

6 

-. ?96“45 

A 

-.632525 

C 

-.584662 

1° 

3 

.6333? 

3 

-.j 67 551 

.22:54 

7 

** 

-. 0 1 3 7 P  0 

4 

.164595 

4 

-.081173 

9 

1. 72246 ” 

c 

* 

-.387:21 

A 

-1.02857F 

C 

-.436574 

19 

2 

1.75868 

7 

-.141261 

•  37 175 

7 

-. C?  ?  s2? 

4 

1. 161351 

4 

-« 3172°7 

A 

-1.6:3411 

0 

-. 33  4  Cl° 

19 

1 

8.48336 

3 

1.0  39 .6  r 

.5f  892 

3 

-.171P04 

A 

-2  •  86 2 

C 

-  »  1?1  ci-  4 

GAK/^ATI -/72-7 


TAFLF  I 


•  rJ?;^Trif.MTC  F0;?  3rsT  -CKHITlOtoAL  ESTIMATION  OF  ThE 
LOCATION  Affl  SCALE  FA”8ME7CF?  OF  THC  CALCHY  CIST9IDLTTCN 
(WITH  ADDITIONAL  C£NC05ING  FROM  AnCVF)~ 


**  LOCATION  ** 

K  H  M°f  I  COcF. 


**  SCALE  *# 

MSE  I  CPFF, 


20  16 


.1256°  3 

‘  A 

5 

6 
7 
A 
9 

1C 
11 
17 
13 
1 A 

15 

16 

‘  17 

10 
A. 


-. ^10562 
-.024  26  8 
-.  .‘28212 
-.nosec  ? 

.041 .52 
.113576 
.135798 
•  231 6C"7 
•231614 

.  135  etc- 
.112161 
.  0  4 1  •*'  t  o 
- .  1  3  3  6  C  ? 
"  •  2  P  2  f  ? 
-.  J24  76  P 
- . 01056? 
.O.'CJCC 


• 12C57  3 

4 

5 

e 

7 

8 
c 

10 

11 

12 

12 

14 

15 

16 

17 

18 
C 


-•0J5n74 
-.  9254r  7 
-.060  f P4 
-.10461: 
-• 14 j  C91 
-.14671 4 
-.112563 
-. 042H7 
. 042229 
.112574 
.145715 
.14,;  569 
.104610 
. «6??84 
.0254:  7 
.005674 
•  C  0  C  0 1  c 


15 


.12624  3 

4 

5 

6 

7 

8 
9 

1-3 

11 

12 

17 

14 

15 

16 
17 

A 


-.91C6LP 
-.724373 
-.023222 
-  ♦  C  3  8  6  ?  4 
.  J41244 
.112691 
.186642 
.222665 

•  722677 
. 1 86666 

•  1 1  ?  7f  5 
.04132° 

-.00P529 
-.325212 
-.  J43  .-21 
.U72525 


.12074  3 

4 

5 

6 
7 
P 
o 

10 

11 

12 

13 

14 

15 

16 
17 

C 


-.  C.J56P2 
-. u?5443 
-.:&2671 
-.104961 

-.140205 

-.146921 
-.112744 
-. 042428 

♦  342  ■*05 

.112691 
.146«72 
.14 j  23  C 
.10  4  P96 
.  0*9907 
.034119 
.032424 


fwgm i>i.;mnwfBnWPW  U.1 


GAN/MA™/72-7 


7APLE  T 

COEFF ICIFNTC  FOR  59FT  CCNf'ITIONAL  ESTIFATICN  C*  IFF 
L  COAT  ION  AN?  CGALE  CAR AN- T£FS  CF  Tpr  CALCHT  rjSTRIPUTTCN 
(WITH  AOOTTIONAL  C 6 NSC RING  FRCF  A9CVF) 


N  H 


2d  14 


20  12 


**  LOCATION  **  **  SCALE  ** 


MSF 

I 

CCSF. 

M«?E 

I 

CCFF. 

1274ft 

7 

-.910711 

,12164 

3 

-.305526 

4 

-  .9  2461  ft. 

4 

-.025635 

5 

-.923590 

5 

-.  061473 

6 

-.398695 

6 

1C  5  76  0 

7 

•041694 

7 

-.141276 

a 

,114867 

e 

-.14  3^9 

9 

.1195P2 

c 

+ 

-.113664 

1 J 

.2T511F 

1C 

-.042*9? 

11 

.235177 

11 

.04245? 

12 

.18579? 

12 

.113290 

IT 

« 1 1 5  2  C  4 

13 

.1476*5 

14 

•  '3421T? 

14 

.140559 

15 

-.  .0*21” 

15 

,195  “24 

16 

-«C°0744 

16 

.  1 C  5  8  0  4 

A 

.  J11265 

C 

.Cl  0  76  9 

1??C  2 

3 

-.010*4:. 

.12437 

3 

-,  03F  r59 

4 

-.02451 C 

4 

- . C26  ?12 

5 

-.028914 

5 

-,  062  822 

6 

-.093747 

6 

-.1^9169 

7 

,  0  4  2  3  C  7 

7 

-.144640 

8 

.116455 

8 

-• 151 59 1 

9 

.191154 

0 

-.11 65r2 

1  J 

.23*747 

10 

-. u44  329 

11 

.275c2? 

11 

, G42  7?9 

12 

,191676 

12 

*115  019 

12 

,  1 1  7  2  ft  6 

13 

.150C53 

14 

.9*4  3  36*. 

34 

.143777 

15 

-.135714 

15 

.236547 

A 

.329121 

r 

.327167 

GAH/MATH/72-3 


T  Anl  E  I 

CCrFFiri»  NTS  FC°  9F5T  CCNC  T  T I  CN'AL  ESTIMATION  CF  THE 
LCCATION  aro  FC  ALE  PARAMETER'S  rc  THE  CAUCHY  CISTPI9UTIC*1' 
(WITM  A  COITION  A  L  C£NSCCING  F»CM  APCVL) 


N  M 


20  12 


2-  11 


20  10 


+*  LOCATION  **  **  SCALE  ** 


Mc£ 

I 

CG  :F . 

MSE 

I 

CCEP. 

13u  19 

3 

-  •  o  1  o  9  7  6 

.12039 

3 

-.7*6760 

4 

-.*>25113 

4 

027617 

5 

-.029140 

c 

-.cesses 

6 

-«Q-'19?42 

6 

-.113671 

7 

.*■42841 

7 

-.152*42 

r. 

♦117/79 

8 

-.1R9E°5 

9 

.193271 

0 

-.123144 

1  j 

♦241.41 

10 

-,04  77r:2 

11 

.241761 

51 

.  04-7 ’77 

12 

.194240 

12 

•  113cp7 

13 

♦1192F7 

13 

.1551*P 

14 

-.075S79 

14 

•  42  3  ?C  7 

A 

.049336 

C 

.  04941D 

lioE? 

3 

019^62 

.14136 

7 

"w 

-.CG*9’1 

4 

32516* 

4 

- .  C  3  .  I  ‘  1 

5 

-.0291F3 

5 

-.^7 1^5° 

6 

-.:?8c97 

6 

-.124  CEO 

7 

.24 3  CPS 

7 

-.166243 

0 

.ii°?£>7 

8 

-.174961 

9 

» 1 9  4  •.  9 1 

c 

-.13599^ 

10 

.242394 

10 

-. C54644 

11 

.24293c 

11 

.043634 

12 

» 1 9  9  o  8  4 

12 

.12514’ 

13 

.0 *»026 

13 

.665 ’6  5 

A 

.067  «24 

C 

•  069  266 

13('68 

-,01*96 3 

• . 16115 

3 

- .  u  u  7  9  0  3 

4 

-.029164 

4 

0^4212 

5 

-.3  2921 9 

5 

-.032133 

6 

-.a?',74r- 

6 

-.141792 

7 

•04732? 

7 

-.130470 

8 

.11  :l2 

8 

-.2J1413 

9 

.19797? 

c 

-.153349 

10 

.241926 

10 

-.'.*7  26*5 

11 

.24272? 

11 

,ii43245 

I? 

.234640 

12 

.010671 

A 

.n5 695  6 

0 

•  C  7  0  ?  4  2 

GA^/w5TU/72-? 


TA  PLE  I 

CCpFFinTcf:T?  F07  RRST  CCNPTT  ICf-'AL  ESTIHATICK  CF  THE 
ICPATICN  A^'O  ^CALf  FAPAMF.TrPS  CF  THE  CAUCHY  CISTF»IRUTICN 
C4TTM  ADDITIONAL  CENSORING  FPCH  ARCVF) 


#  •  •  4 

N 

>  •  •  •  • 

y 

**  LOCATION  ** 

MSfc  I  CO  5  F , 

*¥ 

MSE 

SCALE 

I 

** 

CCEF  • 

2C 

9 

•  lv224 

7 

•  19035 

w 

-.009279 

A 

-•32555* 

A 

-.040654 

5 

-.029745 

5 

-.o 97727 

6 

-. .9*2CR 

e 

-. 169  ufl 

7 

« u  A  2  9  2  A 

7 

-. 227SC 1 

3 

•  1 1  °67  2 

e 

-.?42?95 

q 

.194552 

c 

-.193606 

1  u 

.242:5® 

10 

-. 388  2" 8 

li 

.476^3* 

31 

1. 094773 

A 

. J176C9 

0 

.025246 

2G 

? 

.13»i5 

3 

-.''•11679 

.22660 

3 

-» .li 2A5 

A 

Q :6C?2 

A 

-•  0492?-. 

9 

-«  j21 6C7 

5 

-.  11B527 

6 

-.oipfL? 

f: 

-.205611 

7 

. .  j  4  3  2  f  ° 

7 

-.274745 

? 

. 12 3  5r  9 

6 

-.2987»7 

9 

•23*257 

c 

* 

-.247341 

10 

.71592? 

1C 

1.1~97. 7 

A 

—  •*J37**7^ 

c 

-,u957?6 

21 

7 

.15392 

3 

-.013345 

.26750 

C 

-.0127°? 

A 

-  •  "*  3  r  2E  F 

A 

-. DR?1?R 

5 

-.03513. 

5 

14*235 

6 

-,ii39<-r 

6 

244 106 

7 

•  94412? 

7 

-.322  :?6 

3 

.12032? 

8 

-.259779 

9 

.020211 

0 

* 

.8547^6 

A 

-«ieoAr? 

0 

-.29?6°4 

cl 

6 

.1*59? 

7 

-.015=15 

. 2925C 

3 

-.014765 

A 

:37?.pf; 

A 

-.56 4222 

5 

- .  1 A  5  6  7  A 

5 

-.15517* 

9 

-.02122? 

P 

-•  27'.  992 

7 

.94562" 

7 

-,7896?1 

9 

1.17  A'»'  7 

a 

.27Cf»7 

A 

-.71457* 

D 

-.494  °4 1 

■ . .  jun-ni.ii.  ■  i» 


Gflp/MflTM/7?-7 


TA°LfT  T 

COEFFICIENTS  FC?  9P.ST  CCNDIT IOMAL  ESTIMATION  CF  TM£ 
LCCATTCN  AND  SCALE  FA9AM5 T«"S  OF  THE  CAUCHY  CICTRI°UTICN 
(WTTM  ADDITIONAL  CENSORING  F&CN  A90VE) 


**  LOCATION  ** 

N  M  MSE  I  CO=F. 


**  SCALE  ** 

MSF  I  CCEF  • 


2n 

5 

3 

-.121772 

.29862 

7 

-.015158 

A 

-.391523 

A 

-•  065c72 

9 

-.055599 

5 

-.159556 

6 

-.137460 

£ 

-.278620 

7 

i,17Scll 

7 

-.•279275 

A 

-.532757 

C 

-.598533 

21 

A 

.3917? 

3 

-.024667 

.29963 

7 

-.21514C 

4 

-  •  0  0  A  1  ?  r< 

A 

-„6SK519 

r 

-.111598 

5 

-.159070 

6 

1.2  2r794 

6 

- .  3  v  A  7  C  9 

A 

-.752676 

C 

-.575154 

n 

CM 

3 

3 

-.66794' 

•31345 

7 

-.6156-13 

A 

-.16° 66? 

A 

-.  j&7  8 A 0 

5 

1 .  ?  3  ?  7  r  ? 

c 

-.28394* 

A 

-1.112:98 

c 

-.467397 

21 

2 

1.9413? 

3 

-.182-25 

.36520 

7 

w 

-.617282 

A 

1.1822?° 

A 

-•?342p9 

A 

-1.71067: 

C 

-. 321 e7? 

20 

1 

9.426T1 

7 

i. :i«vir 

.5Cf 71 

7 

V 

-.162640 

A 

-2.025722 

c 

-.16281*0 

f  o 


GAM/MATH/ 72 -3 


APPENDIX  B 
Table  II 


r 


70 


73 


■wf  in,  ngnw 
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T ApL F  II 

COFFFI^IrNTS  FOP  9EST  CCNCTTICNAl  ESTIMATION  CF  THE 
LTCATICN  AMD  SCALE  FAFAMST5FS  CF  THF  CAUCHY  CIST9I5UTT CN 
(WITH  AtlDTTlCNAL  SYMMETRIC  CENSORING ) 


**  location  ** 

N  H  MSE  I  COEF 


SCALE  ** 

MSE  "l  ~  COEF. 


lu 

6 

.32626 

3 

“•3c2'C4 

•'  ♦  2Q771 

•a 

w 

-.1*30112 

4 

.3329F3 

\ 

4 

-.307114 

5 

.475 :F2 

e 

j 

-.212155 

6 

.4  7  5  :62 

6 

.212155 

7 

.33395? 

7 

.302114 

S 

-.  362  if  4 

e 

.103113 

A 

-  •  ‘J  3  3  j  0  0 

c 

3 ,  0  J  j  C  L  o 

10 

A 

.73622 

4 

.2903U 

.72141 

4 

-.483473 

5 

•  4  9  9  ?  C  0 

5 

-,229c69 

6 

. 49°  7C " 

6 

•  2  2  9  CP  9 

7 

.  1  j  3  c  C r’ 

7 

.437478 

A 

- .  c  3  o  •:  c- : 

C 

.  o  c.  o  •:  ••  o 

10 

2 

.33622 

c 

-» 

.  5  C  3  •:  ?  C 

.51 9 26 

5 

-1. 23?  071 

6 

.F33.tr 

6 

1.233171 

& 

j.  j  3  3  .  t;  C 

r 

.COCCv.0 

11 

7 

.23197 

3 

-.09212.0 

.25747 

i 

-.06761 3 

4 

.00865c 

4 

-.231556 

5 

,299622 

5 

-.275726 

6 

.496199 

6 

—  •(>(!. (tO 

7 

«  2  9  5  e  ?  2 

7 

.275726 

3 

•  0  3  8  6  5  9 

8 

. 231 5n6 

9 

-.05228: 

q 

.067613 

A 

-.o-jcrcc 

r 

-.  JiiCOC  0 

11 

5 

.2<*f47 

4 

-.  066  13  f; 

.27193 

4 

-.345120 

5 

.303367 

5 

-, 289 1°9 

6 

•515459 

6 

.  C  r  j  f  C  Q 

7 

.303387 

7 

.289199 

5 

-.366116 

5 

•345120 

A 

.000 jL" 

p 

.003000 

11 

3 

.2  <i2Q? 

5 

.2365F4 

.26217 

5 

-.932166 

6 

.52*r92 

6 

.OOOOCO 

7 

.226954 

7 

.932166 

A 

-  •  'j  0  0  c  c  c 

C 

.nc  octet 

7? 
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TA9LE.  II 

COEFF  ICIf. NTS  FC9  EFST  CCNPITICNAl  ESTIMATION  C*  TFC 
LOCATION  AMO  SCALE  PARAMETERS  OF  THE  CAUCHY  CISTRIDUTICN 
(WITH  AT9ITICNAL  SYM*FTRIf  CENTRING) 


**  L03ATI0U  **  **  SCALE  ** 


N 

M 

msc 

I 

CCEF  • 

MSF 

I 

CCEF. 

11 

7 

. 22673 

A 

-.085or4 

•  21C95 

4 

-.184503 

5 

.084779 

5 

-.237471 

6 

.?9  8  Of  8 

6 

-»  2056G2 

7 

.412411 

7 

.-j  core  a 

3 

.293868 

8 

•  205  60  2 

9 

•  n  8  4  7  7  9 

o 

.237471 

10 

-  •  039854 

10 

.1845C ? 

A 

-.OgQCfC 

r 

-.000:10 

13 

5 

.23261 

5 

-.9*4512 

•  2  2  P  4  rJ 

C 

✓ 

-.524726 

6 

.310761 

6 

-.22°124 

7 

.427501 

7 

*  C  0  0  p  2  u 

8 

.310761 

8 

.  229124 

0 

-.024912 

o 

.524726 

A 

•  0  j  0  C  C  v 

*3 

.00  0  CO? 

13 

7 

.23270 

6 

.  2  3  5  4  r  5 

.  74462 

e 

-1.17? 245 

7 

.4?oiE3 

7 

•  Ct.  0wL4 

8 

•  2  3  9  4  f  « 

8 

1.172245 

A 

-.COD  it  0 

(? 

•  j  C  G  C  £  0 

14 

ir 

. 199? j 

7 

-.02917*' 

.18722 

3 

-.025067 

4 

-.34n524 

4 

-.lt0iB7 

5 

•03171? 

5 

-.19450  2 

6 

.195:ir 

F 

-.215877 

7 

.342931 

7 

-.0«6717 

8 

.342931 

8 

.096717 

9 

.199011 

9 

•  215  fl77 

10 

.031712 

10 

.194902 

11 

-  •  9  4  J  5  2  4 

11 

.  1  j  0 1 P.  3 

12 

-.3291?'' 

12 

.o?5 re  7 

A 

-.0  30  H-G 

r> 

W> 

-.ccorco 
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7A*LF  II 

CPlFFICI^MTS  FOR  5EST  CGNDIT I C*JAL  ESTIKATICN  C F  THE 
LCCATTCN  AMO  SCALE  FAC AMETrFS  OF  THE  C AL'CHY  CI'TRIPUTICM 
(WITH  ADDITIONAL  SY^MFTRIC  CENSORING) 


N  M 


15  11 


15  o 


15  7 


15  5 


**  LOCATION  **  **  SCALE  ** 


M?t 

I 

CCee  * 

KSE 

I 

CCEF. 

18231 

3 

-.02A13f> 

.17153 

1 

w 

-.018981 

tt 

-.029711 

9 

-.077236 

5 

.3327*6 

5 

-.160*27 

6 

.123*31 

6 

-» 2C6229 

7 

•  259219 

7 

- .  153  C95 

5 

.  3  2  5  EE  ^ 

8 

-.030000 

<3 

.259319 

9 

. 15 o  096 

<  0 

.  122  331 

10 

.296229 

11 

.3027*6 

11 

.160*23 

12 

-.320711 

12 

.077336 

13 

-.02*130 

12 

.C189~l 

A 

.0CC-3U 

C 

. GOD  0l 0 

ir.y?P 

t* 

-.u7e;°2 

•  1 73F0 

* 

-. 106028 

? 

.••’33773 

5 

-.161909 

6 

.x2EPf o 

6 

-.20 3221 

7 

•  2  7  *  7  f  C- 

7 

-.151670 

5 

.3*25* E 

8 

.  0  C  j  t  C  0 

9 

.27*780 

9 

.161670 

10 

.126552 

1C 

.239221 

11 

.033273 

11 

.161909 

.12 

- .  0  7  6  C  9  2 

12 

« 1  u  6  0  2  9 

A 

-.3.33  :co 

C 

•  QC  OCiiO 

,19pe5 

c 

«<* 

-.09**56 

. 18*29 

5 

-.313609 

6 

.  1 32  31  7 

6 

2195*0 

7 

.?s*  °re 

7 

-, 16n 12* 

9 

.35*6*1 

8 

.000303 

9 

.  ?e*4  e**  ? 

0 

,  1 6  0  1 2  A 

10 

.172217 

10 

.2195*3 

11 

-.09**96 

11 

,  212  6C  9 

* 

•n 

•  C  *  u  ij  0  Cl 

c 

-.DCO-CO 

6 

« 0 26  7?6 

.22037 

6 

-«693cll 

7 

.291  98  * 

7 

-.188C73 

3 

. 2627P1 

8 

•  jO.UtC 

9 

,291sr* 

c 

* 

.188-73 

13 

.  'J 7  6  ?  2  6 

10 

.693911 

A 

-•jCS--CC 

r 

*  #  0  C  j  'CO 
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TAPLE  T I 

00£ FPTCIFNT?  FOR  B^ST  CCNTITICNAL  ESTIMATION  CF  Thc 
LOCATION  AMO  SCALE  PARAMETERS  C*  THE  CAUCHY  CISTRI°UTICN 
(HITH  ADOITICNAL  SYMMETRIC  CENSORING) 


N  M 


16  12 


16  JO 


**  LOCATION  **  **  SCALE  ** 


MSE 

I 

COSF. 

M4'E 

I 

COcF . 

19351 

7 

.219-fF 

.33635 

7 

-1. 4 J4o20 

8 

.361 cRG 

8 

—  .00CCLC- 

9 

.319JC5 

c 

1. 494620 

A 

-.000  etc 

C 

.  ooc  :ro 

16709 

3 

-.02014? 

.15819 

3 

-.014624 

A 

-.C370F7 

4 

-.00047? 

5 

-.013592 

c 

•* 

-.131767 

6 

.07266-2 

6 

-.136477 

7 

,2i  0627 

7 

-.1  221Q 

8 

.296514 

e 

-.070793 

a 

* 

.296514 

c 

•* 

.073798 

n 

•  2C  3  c2 7 

Id 

.172219 

n 

.  773062 

11 

.18*477 

12 

-.Gl'iSo? 

12 

.131267 

13 

-.027:67 

13 

. 06347? 

14 

-.•J20142 

14 

.014524 

A 

,933000 

n 

-•uOCjtC 

16  °67 

A 

-.357479 

. 15952 

4 

-.0822*3 

5 

-•0124C3 

5 

-.122222 

6 

.  375214 

6 

-.187634 

7 

.234147 

7 

-.172627 

3 

.331521 

8 

-.071245 

9 

.3215?! 

c 

.0-71245 

10 

.234147 

10 

.173627 

11 

.075214 

11 

. 1 87  ?74 

12 

-. 0134( ’ 

12 

.132222 

l7 

-.067470 

13 

.082362 

A 

.  3  3  0  36.  C. 

C 

•  0  C  j  C  P  0 
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T AnLE  II 

COEFFICIENTS  FOR  £rST  OCNPIT  TONAL  fSTIMATICN  C*1-  TH^ 
LCOATICN  Af'3  SCALE  PARAMETOFS  OF  TPF  CACC^Y  ricTRI5'JTICK 
<WITH  ACOITICML  SYMMETRIC  CENSORING ) 


**  LOCATICN  **  **  JC^LE  *+ 


N 

M 

'  MFc 

I 

CC5F. 

KEF 

I 

CCLF. 

16 

fi 

.1744C 

5 

-.1  31441 

• 16670 

c 

* 

-.247265 

6 

.276711 

6 

-»195rc’3 

7 

. 211 2f  2 

7 

-.160517 

6 

.’HA?? 

e 

-. J74253 

q 

.211423 

0 

. C742r  7 

10 

,211262 

lu 

,  1 A  r  r 7 

11 

.0 7*711 

11 

,  1 9  5  C  f  3 

12 

-.101441 

12 

.247235 

A 

.00 011  c 

C 

-  •  J  J  J  C  l  o 

16 

o 

.1776? 

b 

-.M3:  F." 

.16991 

e 

-•  543  ”45 

7 

. 2 1 7  ci C 

7 

-,2J2c5f. 

3 

.222144 

6 

2 

a 

.120144 

C 

•» 

.  0«.?Pi  •» 

10 

.217?jf 

ID 

.2*-2  66j 

11 

-.073C6C 

11 

.  54  9  54  5 

A 

•  f  ;  ?  cc  r 

C 

0.0  >•  9  : 1  5 

16 

A 

.17612 

7 

» 17?264 

•  25614 

7 

-1. 071^3 

6 

.32  IP  20 

e 

-,11-Of.O 

0 

•32162F 

0 

.  1 1 :  c  l  ;• 

10 

.176764 

10 

1.  071162 

A 

•0JCCC3 

C 

Cl.uOG*0  0 

15 

2 

.1*145 

6 

.  5  J  o :  C  0 

,49464 

6 

-2.226576 

q 

•  50031  'j 

c 

2.726576 

A 

j.oc c:Ln 

C 

.  0  0  3  C .  0 
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TfiHF  II 

CCrFFICTLNTS  FOR  SFST  CCNFITTCNAL  5STIHA7ICN  OF  THr 
ICCATTCN  A Nn  SCALF  FA^AmhI'IRS  OF  THF  CAICHY  r IS! PI°L'TT  CK' 
(WITH  ADDITIONAL  S YNM^TRIC  CENSORING) 


ICATIOM  ** 

I  C05F. 

** 

MSE 

SCAL 

I 

r  ** 

OOFF. 

? 

-•315c48 

.14677 

■2 

-.-11343 

4 

-.037775 

4 

-.C47P78 

5 

-.322251 

5 

-•i*»7G5i 

6 

.03527? 

6 

-.163519 

7 

.142717 

7 

-.175151 

3 

.249  75/- 

8 

-.11501* 

g 

.233431 

9 

.  0  0  ?  0  0  9 

1J 

.  2*«F  754 

1  0 

.115115 

11 

.143717 

11 

.175151 

12 

.039273 

12 

.163919 

13 

-.022251 

13 

.137651. 

14 

-.322775 

14 

,  °47  f  7  8 

15 

-.11 6 5* 8 

15 

. J11243 

A 

.  j  c  o ;  c  r 

C 

. JGGCC0 

4 

-.05927? 

.14769 

4 

-.-64C' n 

n 

-.022774 

5 

-.1-3151 

ft 

.  3  4  0 1 P  9 

6 

-.164791 

7 

.  145  925 

7 

-.176132 

3 

.249222 

8 

-.115673 

9 

.292*20 

o 

* 

-.OCOCwC 

10 

.24???? 

10 

.116673 

11 

« 1 4  6  °  2  5 

11 

.17613? 

12 

.040189 

12 

.104791 

1  ? 

-.02227/. 

12 

.103161 

14 

-.03937? 

14 

. 064 er 3 

A 

*  C  C  j  j  t  c 

C 

-  •  j  r  o  c  -  o 

9 

-.103543 

.15?6C 

c 

✓ 

-.197322 

6 

.34225" 

6 

-• 16°417 

7 

.15075  3 

7 

-.181 354 

3 

.255278 

8 

-.119172 

o 

.3  Cl  712 

c 

•  #  3  j  j  i.  c 

10 

.255.878 

10 

.11917? 

11 

.15978 2 

11 

.181364 

1? 

•  0  422*0 

12 

,  16941  7 

12 

-,liOcM 

12 

.197332 

A 

- .  o :  o  •- 1 r 

C 

-.jPlttO 

N  M 


M  «*£ 


17  13 


.15440 


17  il 


,1!  646 


17  9 


.16055 


79 
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TA°LF  TI 

COEFFICIENT*  FOR  CC*T  *CND JTI CK'AL  ESTIMATION  CF  TH* 
LCCATJCN  A^n  SCALE  FA°AMET5FS  CF  THE  CAUCHY  DISTRIBUTION 
(WITH  ADDITIONAL  SYMMETRIC  CENSORING) 


♦  •  • 

N 

V 

*  * 

MSE 

LCCATION  ** 

I  CO  =  F. 

*¥ 

MSE 

SCALE  ** 

I  COEF. 

17 

7 

.1642° 

b 

-.••'■76144 

.16819 

6 

-.44°  794 

7 

.136  3  (-5 

7 

-.15774? 

s 

.26433  5 

8 

-.13019° 

9 

.210326 

9 

,  v  J  'J  r  0  0 

13 

,  264  °3  5 

1C 

.13019° 

11 

.156  US 

11 

.197742 

3  2 

-.376,14  4 

12 

. 440794 

A 

- .  0  r.  C  C  C 

r 

— «  C  u  J  C  0  3 

17 

5 

•  1 6  5  ?  8 

7 

.074965 

.20976 

7 

-.847475 

3 

,26" 37 3 

8 

-,159c 9 7 

o 

.212267 

q 

.  i f  C .  c 

10 

. 2C  9  ,72 

10 

.15 9293 

3.1 

•  J  7  4  5  6  5 

11 

. 947475 

A 

.  0  3  u  j  L  r 

C 

.vJ'JC,0 

17 

3 

.1(527 

3 

.  347  4  rr 

.  3  ?  r  C  6 

8 

-1. 63?4?7 

9 

,713  '27 

Q 

-  •  j  G  T  ■'  >-  li 

10 

,34346? 

10 

1.63?4«7 

A 

-.C*l  j'.tC 

c 

•  0  t*  -  .‘t  » 

1*) 

14 

.1  434° 

3 

-.01437? 

.13687 

3 

-. :r8E°7 

4 

-.270790 

4 

-. u38257 

- .  0  ?  6  f  0  5 

5 

-.088544 

5 

.016373 

6 

-. 142010 

7 

.099335 

7 

-.167494 

3 

,195378 

ft 

-. 1?8745 

9 

•  2  6  •*'  2  3  4 

S 

-.052664 

U 

«  2  G  n  3  7  4 

10 

.  3*3884 

11 

,19533" 

11 

•  1 3  "  3  4  5 

12 

,  0  9  9  2  3  K 

12 

.167494 

1? 

•'16272 

33 

,  1  h  2  C 1 0 

14 

•*#  y ? 6 6 ^ ^ 

14 

.  08  8  54  4 

15 

-, 02:399 

35 

.G78257 

1 6 

j 1 4  77? 

16 

.  JO 5 9° 7 

A 

.0  2CC(  C 

C 

-.tCGCCO 

?r 
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TABLE  II 

CC£FPTCTt  MT°  FC9  BEST  CCNTITTCMAL  r<;TIMATICN  OF  THE 
LCCATTCN  Af"  ^CALc  F  A  *  A  M  E  T  c  F.  S  OF  THE  CAt.THY  CISTRI°UTICN 
(HIT w  AODITm.AL  ? Y ^ ‘*1  u  T?IC  CENSORING) 


N  H 


16  12 


13  lb 


19  6 


LOCATION  **  **  SCALE  ** 


;  C-- 

I 

COE*. 

MSE 

I 

CCEF. 

1451t« 

4 

-.052098 

. 13752 

4 

-.051812 

5 

- . 02 5  ?C « 

p 

+ 

-.088841 

6 

•  F 16  77? 

6 

-.  142c7  0 

7 

.lOOFCr, 

7 

-.16*198 

6 

,197798 

a 

-.138=4? 

0 

.283545 

0 

-.05411? 

1  0 

« 267545 

lu 

.054119 

11 

,127  79-? 

11 

.138=4? 

12 

.12099c 

12 

.1^9393 

13 

..*16  772 

1? 

.142570 

14 

-.028/(0 

14 

.0*3*41 

15 

-. J32  :?* 

35 

.  051*32 

A 

. f :  3 :  c  e 

C 

.1 
r  . 

o 

-> 

•  *i 

• 

1 

14*64 

5 

-. -75757 

•  1 4  C  9  7 

5 

-.159215 

6 

,ni7c2? 

6 

-.145533 

7 

« 1 97  8ct 

7 

-.171954 

A 

.2033?? 

a 

-.14911 C 

9 

.270-027 

c 

* 

-.055273 

13 

V  . 2?0c?7 

10 

.-55273 

11 

« 2  0  7  *  C  ? 

11 

.14231 0 

12 

•  1  u  7  8  5  6 

12 

.  l“'l  =  24 

13 

,017222 

4  7 

•  w 

.145538 

14 

-.  )95?C7 

14 

. 15°  21 5 

A 

.COCObC 

C 

.  L'jOCt  C 

1 15235 

6 

- « 096677 

.15179 

6 

-.35*577 

7 

.1079^ 

7 

-.  1.8  3  524 

8 

.  2  19  873 

8 

-.  151984 

9 

•  2  7  fl  9  G  G 

c 

-.059223 

10 

, 27  8  90  v 

10 

. ^59273 

11 

.2  05  B71 

31 

.151=84 

12 

,ir-72ff 

32 

.18:524 

13 

-.8=9677 

33 

•  ?c  *  877 

A 

•uB5wL' 

D 

-  .  u  l  J  7  'j  0 

**  5 
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TflrLF  xi 

CCrFFICIfK^?  FOR  5FST  CCNrTTICNAL  rSTTf-*AT ICN  CF  TFr 
LCCATTOK  AMO  SCAL£  PARAM“TrF.S  CF  ThF  CALC^Y  CI5TRI°L,TTCN 
(WITH  aooiticnal  sywm-tric  CEN*C9!NC) 


**  LOCATION  **  SCALF  ** 


N 

W 

M5£ 

I 

ccef. 

M^F 

I 

CCFF. 

IS 

6 

.IFA-jF 

7 

.002:92 

.17941 

7 

-.F*8F,S9 

a 

•217c 7P 

8 

-.  176919 

9 

.2*2?7C 

c 

-• 06° lu  1 

10 

« 2  8  7  77" 

10 

.06°  1-1 

11 

» 2 12  87R 

11 

.176919 

12 

. C  '729? 

12 

•68*269 

A 

/:?ai 

C 

i  u  t  'j  r  i>  *i 

IS 

A 

.1 5 A' F 

a 

.21C717 

.  ?9-r  C3 

0 

-i.249C° j 

9 

.  2  0  3  0  6  ? 

c 

-.094011 

u 

.  2  *3  ?  €  e  ? 

1C 

.c°4:n 

u 

.21021'’ 

11 

1. ?4°  r9  0 

A 

.  c  0  v  2  L ; 

0 

.  o  j  :*  c  w  3 

id 

7 

.If  709 

a 

.  5  ;  j  r ; 

.  4  c  2:  /-  9 

c 

- 2. 66fl  4^  A 

*  .* 

X  •• 

.  5 : :  *  r « 

1  J 

2.66  4r  4 

A 

: ,  ?  c  o  j  c  r- 

C 

. » :■  •:  :*c 

19 

15 

•? 

-.31?’76 

.  1PP21 

*1 

V 

-.007227 

H 

- .  :•  2  7 1 9  4 

4 

-. .31769 

5 

- .  j  2  ?  1  f  « 

5 

-. .732" 8 

6 

.  0  J12CP 

6 

-.12216? 

7 

.965^  7 

7 

-.154747 

6 

,  i5j5°r 

8 

-•147151 

9 

.227f°? 

c 

-.09,45° 

i; 

.252097 

10 

;  jOvCCO 

11 

.222092 

11 

. T9? 469 

12 

.15'9°' 

12 

•  147 !rl 

13 

•  0  0  5  0  4  3 

12 

.154747 

14 

« «.  3 1 2  (  C 

14 

.122162 

15 

-•P28ir«* 

15 

.  n7  3?r s 

16 

-  •  *»  ?7  1  C(, 

16 

.271659 

17 

-.012270 

17 

.  ,'.r?227 

A 

-.riMr 

C 

.  c  c  j  •:  w  o 
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T  A  5LE  II 

CCLFPICIFNTS  FO?  rE?7  CC^.TTTTONU  F ST  I Y AT I C  N  CF  TK 
LCCATICfi  AHD  SCALE  PA^AMET  EPS  3F  THE  CAUCHY  C  I^TRI?IJTI  CM 
(WITH  ADOITKMl  T9IC  C  EhSCFING ) 


N  H 


19  12 


19  11 


i°  g 


**  LOCATION  +* 
msc  i  no?c. 

** 

msf 

SCALE 

I 

CCEF. 

-.o.»572? 

•12P68 

4 

-.041661 

K 

-.728216 

c 

-.Q733«7 

6 

. ? "1 27  A 

6 

-.122525 

7 

.766656 

7 

-.155.242 

<5 

.15221  1 

9 

-♦147f77 

g 

.226  -PK 

c 

-.  Aqi’62 

17 

,255'!A? 

1C 

•  o  C  C  v  G 

11 

.226  Cf  5 

11 

. 79 3 752 

12 

.152251 

12 

.14? *27 

12 

.766656 

13 

.155242 

14 

.071274 

14 

.122525 

IF 

-.'•2®  216 

15 

.373*,»7 

16 

-.  1-45  722 

16 

.041 p5 1 

A 

.00*765 

r 

. ODOCuO 

l^a??  5 

. 171  ir 

5 

-.125777 

6 

.ooigr 6 

6 

-.124442 

7 

. 76967? 

7 

-.157665 

e 

.156134 

e 

-.150217 

q 

.21156 ? 

Q 

-.092,,f  3 

12 

.?51 431 

15 

•»  -  r. r  r  •  n 
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T  A  PIE  II 

CCCFFIOTffcTS  PC'5  EEPT  5CKDTTICNCI  FSTIMATICN  CF  TWg 
LCCATIOP  PCALC  FA9AM-T9RS  QF  TFE  CALOY  CISTRI’UTICN’ 
(t-'ITM  A’OITTCNAL  SY^TRIC  CENSGCINF) 


N  P 
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**  LOCATION  **  **  FCALF  ** 
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TA51.E  II 

CC  7  pF 10 1 rNT F  FC9  ’’EFT  CCNPIT  IONAL  FSTIKATICN  CF  THr 
LCCATTON  M‘D  SOALE  cAnAKETERS  0*  t>£  CALOY  CIS7RI°UTIC" 
(WITH  AOniTlONAL  S^METRIC  OEK^OPING) 
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APPENDIX  C 
Computer  Program 
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PROGRAM  C03FFS  (INPUT,  OUTPUT) 

C  THIS  PROGRAM  COMPUTES  AND  TABLES  THE  COEFFICIENTS  FOR  BEST 
C  LINEAR  ESTIMATION  OF  THE  LOCATION  AND  SCALE  PARAMETERS  OF 
C  THE  CAUCHY  DISTRIBUTION  FOR  SAMPLE  SIZES  OF  5(1)20  tfTM 
C  ADDITIONAL  CENSORING  FROM  ABOVE.  MUST  READ  IN  THE  EXPECTED 

C  VALUES  AND  COVARIANCES  OF  THE  ORDER  STATISTICS  FORMAT  (6F12.6). 

.  DIMENSION  EXP(20,20),C0V(20,20),A(20,20),3(20),X(20),RL(20), 

1SE(20,20) ,AS(20,20) ,BS(20) ,XS(20) ,RKC(20) ,SS5(20) 

C  NNsJ'lAXIMUlI  SAMPLE  SIZE  Mf^MINIMUM  SAMPLE  SIZE 
NN=20 
MM=5 

DO  5  N=MM,  NN 
IJIssN-2 
NRrsli/2 
NM=(N+l)/2 

HEAD  101,(BXP(N,I),I=3,NM) 

DO  4  1=3, NR 

4  EXP(N,N+1-I)=EXP(N,I) 

5  CONTINUE 

DO  100  N=KM, NN 

NIsN-2 

DO  8  1=3, NI 

8  READ  10l,(C0V(I, J), J=I,NI) 

DO  10  1=3, HI 
DO  10  J=I,NI 
10  COV(J,I)=COV(I, J) 

-  C  FILL  THE  A  MATRIX 
A(1 ,1)=0.0 
A(l,2)=0,0 
A(2,l)=0.0 
A(2,2)=l,0 
DO  15  1=3, HI 
A(I,1)=1.0 

15  A(I,2)=-2XP(N,I) 

DO  16  J=3,NI 
A(I, J)=1,0 

16  A(2,J)=-EXP(N, J) 

DO  17  1=3, HI 

DO  17  J=3,NI 

17  A(I, J)=C07(I, J)+EXP(N,I)¥EXP(N, J) 

C  FILL  Ti!S  B  MATRIX 

B(l)=l.o  : 

DO  18  1=2, MI 

18  3(l)=0,0 

C  COMPUTE  THE  COEFFICIENTS  OR  CENSOR  AND  COMPUTE  THE  COEFFICIENTS 
RifeN-4 

DO  30  ICEN=1,N4 

M=N-3-ICEN 

NP=Mf2 

CALL  KTX£L(NP,A,B,X) 

DO  20  1=3, NP 

20  RL(I)=X(I)  ’ 

RK=X(2) 

CALL  SSL( N , M, NP, RX, RL, KXP, COV, S3) 
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FILL  THE  AS  MATRIX 
DO  25  1=1 fN4 
DO  25  J»1,N4 

251mJL°CSx”,+2>tEXP(S'I+2,*E;tf'(!''J+2) 

DO  26  Isl,K4 

26  BS(I)sEXP(NfI+2) 

THE  COEFFICIENTS  FOE  THE  SCALE  PARAMSTH) 

CALL  KrXEL(»jr,AS,BSfXS) 

DO  2?  1=1, NT 

27  RKC(I+2)=XS(I) 

C^SEL(N,MfNP,ltOfRKC,EXP,COV,SES) 

DO  28  10=3, up 

28  SUI®=SUI-ID«i!lKC(ID) 

RDfcSUKD 

100  cMriMCTUN,M'SE,aL,SES,R:'tc,ili'* 3lllD> 

101  FORj'IAT(  6?12«  6) 

STOP 

END 
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SUBROUTINE  PRINT1  ( N, M, SSL,  COFL.SBS ,  COFS , A ,  13,  D) 
DIMENSION  SEL( 20,20), C0FL(20) , SES(20,20),C0FS ( 20) 

THIS  ROUTINE  PRINTS  THE  COEFFICIENTS  AND  MSB  FOR  SINGLE 
CENSORING  FROM  ABOVE. 

N=S AMPLE  SIZE 

SAMPLE  SIZE  AFTER  CENSORING 
SELsMSE  OF  LOCATION  ESTIMATOR 
COFLs  COEFFICIENTS  FOR  THE  LOCATION  ESTIMATE 
SES  =  KSE  OF  THE  SCALE  ESTIMATOR 
COFSs  COEFFICIENTS  FOR  SCALE  ESTIMATE 
As  CONSTANT  FROM  LOCATION  ESTIMATE 
Da  CONSTANT  FROM  SCALE  ESTIMATE 
IF(N.GT.5)G0  TO  5 
IPAGEs39 
PRINT  20 

PRINT  21 ,N,M,SEL(N,M) ,13, C0FL(3) ,SBS(N,M) ,I3»COFS(3) 
M?«M+2 

IF(M,EQ.1)G0  TO  4 

PRINT  22, (I,COFL(I) ,IC0FS(l),Is4,M2) 

4  PRINT  23, A, D 
K«37 
RETURN 

5  RKsM 

RKsK 

RLINEsRI-^2.0 
IF(RLINE.GE.)GO  TO  10 

8  PRINT  21,NlM,SSL(N,M),I3,COFL(3),SSS(N,M)I3»COr'S(3) 
M2=M+2 

IF(M.EQ.l)go  to  9 

PRINT  22, (I, CO?L(I) ,I,C0FS(I),I=4,M2) 

9  PRINT  23, A, D 

KsX-M2 

GO  TO  16 

10  ISKIP=K+3 

DO  11  ISsl.ISKIP 

11  PRINT  24 
PRINT  25.IPAGS 
IPAGEs  IPAGSfl 
PRINT  20 

K=41 

15  GO  TO  8 

16  IF(N.E3.20)GO  TO  17 
RETURN 

17  IF(M.EQ.1)G0  TO  18 
RETURN 

18  isnp=K+3 

DO  19  I3sl,ISXIP 

19  PRINT  24 
PRINT  25,IPAGS 
RETURN 
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20  F0RMAT(1H1 , 1 5X, i 3HGAK/KATII/72-3, ///,4lX,?HTABLE  I, //,19X, 5HC0SFFI 
1CIENTS  FOR  BEST  CONDITIONAL  ESTIMATION  OF  THE,/,17X,5&lLOCATION  AN 
2D  SCALE  PARAMETERS  OF  THE  CAUCHY  DISTRIBUTION , / , 25X ,  39i!( WITH  ADDIT 
3I0NAL  CENSORING  FROM  ABOVE) ,//, 1 5X , 59( *.*),/*  28X ,  14H**L0CATICN** 

4, 16X , 11H**  SCALE  ** , / , 1 6X , 1HN , 3X,1HM, 5X , 3HKSE, 5X ,1HI,4X , 5HC0EF. , 10 
5X, 3HMSE, 5X, 1HI, 4X , 5HC0SF . , / , 1 5X , 59(* . * ) ) 

21  F0RMAr(lH0,l4x,I2fl4tFl0.5,l4,F11.6lFl3.5,l4lF11.6) 

22  FORMAT(>3XlI2,F11.6#15X,:2tF11.6) 

23  FORMAT ( 34X , 1 HA , FI 1 . 6 , 1 6X , 1HD, FI 1 . 6) 

24  F0RMAT(IH  ) 

25  F0RMAT(42X,I3) 

END 


SUBROUTINE  SEL(I.\M, UP, RK,RL, EXP, COV,SB) 

DIMENSION  RL(20) ,EXP(20,20) , COV(20,20) ,S3(20,20) 

C  THIS  ROUTINE  COMPUTES  THE  MSS  FOR  THE  LOCATION  AND  SCALE 
C  ESTIMATES 
SUMluO.O 
SUM>0.0 
DO  10  1=3, NP 
SUM2=0, 0 
DO  5  J=3,X? 

5  SUM2=RL( I ) *RL( J ) *COV(I ,  J )+SUM2 
SUM1=SUM1+SUM2 

10  SUM3=-RL(l)*EXP(N,I)+SUM3 
FA=SUM3-RX 
FA2=FA*FA 
SE(N,M)*SUMi+FA2 
RETURN 
END 


SUBROUTINE  MTXEL(NP,A,B,X) 

DIMENSION  A ( 20 , 20) , B( 20) , X( 20) , PIV( 20) , C( 38 , 38) 

C  THIS  ROUTINE  IS  A  MODIFIED  VERSION  OF  tffXEQ-KATRIX  EQUATION 
C  SOLVER  SUBROUTINE,  COMPUTER  SCIENCE  CENTER,  WRIGHT-PATTER30N 
C  AFB,  OHIO 

C  TO  SOLVE  THE  LINEAR  SYSTEM  AX=B 

DO  10  J=1,NP 
DO  10  1=1, NP 
10  C(I,J)=A(I,J) 

NPJsNP+1 
DO  20  1=1, NP 
20  C(I,NPJ)=B(I) 

NP1=NP+1 
NPKeNP+1 
DO  120  1*1 f BP 
IP1:J+1 
ATPE=0, 0 
DO  40  J=I,NP 

IF  (ABS(C(J,I))-ATPE)  40,30,30 
30  ATPE=AB3(C(J,I)) 

IPIV«J 
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40  CONTINUE 

IF  (ATPE)  210,210,50 
5C  DO  60  JsIPl.NPK 
60  PIV(J)=C(IPIV,J)/C(IPIV,I) 

IFROMsNP 

ITOsrNP 

70  IF  (IFROM-IPIV)  80,100,80 
80  RM=-C(IFR0I1,I) 

DO  90  J=IP1 , NPK 

90  C(IT0,  J)=C(IFROM,  J)+RM*PIV( J) 

IT0=IT0-1 
100  IFR0M=IFR0M-1 

IF  (IFROM-I)  110,70,70 
110  DO  120  JsrIPijNPX 
120  C(I,J)=PIV(J) 

I=NP 
130  IP1=I 
1=1-1 

IF  (I)  160,160,140 
140  DO  150  JsNPl.NPX 
DO  150  LsIPl , NP 

150  C(I, J)=sC(I,  J)-C(I,L)*C(L, J) 

GO  TO  130 
160  NPJ=NP+1 

DO  170  1=1, UP 
170  X(l)=C(I,!!PJ) 

180  RETURN 
210  PRINT  1001 

1001  FORKAT(  37HODET(A)=0  IN  CALL  TO  SUBROUTINE  KTXSL) 
RETURN 
END 
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PROGRAM  COEFFD(  II! PUT ,  OUTPUT  ) 

THIS  PROGRAM  COMPUTES  THE  COEFFICIENTS  FOR  THE  CONDITIONAL 
BEST  LINEAR  INVARIANT  ESTIMATION  OF  THE  LOCATION  AND  SCALE 
PARAMETERS  OF  THE  CAUCHY  DISTRIBUTION  FOR  SAMPLE  SIZES  5(1)20 
WITH  ADDITIONAL  SYMMETRIC  CENSORING.  MUST  READ  IN  THE  EXPECTED 
VALUES  AND  COVARIANCES  OF  THE  ORDER  STATISTICS,  F0RMAT(6F12,6) 
NNsMAXIMUM  SAMPLE  SIZE  KMsMINIMUM  SAMPLE  SIZE 
NN=20 

DO  5  N*KM,  NN 

NIaN-2 

NRsN/2 
NJfc(N+l)/2 

READ  lOl.CSXPCN.D.I^.JJM) 

DO  4  1*3, NR 

4  EXP(N,N+1-I)=-EXP(H,I) 

5  CONTINUE 
DO  100  N=MM,  NN 
NIcN-2 
DO  8  1=3, NI 

8  READ  101,(C0V(I,J),J=I,NI) 

DO  10  1=3, NI 
DO  10  J=I,NI 
10  COV(J,I)=COV(I,J) 

C  FILL  THE  A  MATRIX 
A(1,1)=0,0 
A(l,2)=0,0 
A(2,l)=0.0 
A(2,2)=l,0 
DO  15  1=3, HI 
A(I,l)=1.0 

15  A(I,2)=-EXP(N,I) 

DO  16  Jsr3.NI 
A(1,J)=1.0 

16  A(2,J)sr-EXP(N,  J) 

DO  17  1*3, NI 

DO  17  J=3,NI 

17  A(I,J)=COV(IfJ)+EXP(N,I)*EXP(N,J) 

C  FILL  THE  B  MATRIX 

B(1)=1,0 
DO  18  1=2 ,NI 

18  B(I)=0.0 

C  COMPUTE  THE  COEFFICIENTS  FOR  THE  BASIC  CENSORED  SAMPLE 
KsN-4 
NP=Mf2 

CALL  MTXBL(NP,A,B,X) 

DO  19  1=3, UP 

19  RL(I)*X(I) 

RKssX(2) 

CAI.L  SSLD(K, M.RK ,  RL, EXP,  COV, SSD) 

C  FILL  THE  AS  MATRIX 
R4-N-4 
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DO  38  1=1 » II 
38  RKC(I+2+ICEN)kXS(I) 

CALL  seld(ii,m,i.o,r:<c,exp,cov,sbs) 

lB=3fics:i 

ITsIB+M-l 
SUMDbO.O 
DO  41  XDsID.IT 
41  SUKD=SU!1D+RKC(ID) 

RKDsSUMD 

50  CALL  PRINT ( N ,  M,  SED,  RL ,  SES , RXC , RX, 18 ,  RXD) 

100  continue; 

101  F0RHAT(6F12.6) 

STOP 

END 


SUBROUTINE  S3LD(  N ,  H.  RK ,  RL ,  SXP ,  COV,  SED) 

DIMENSION  RL(20) , EXP ( 20 , 20 ) , CO V( 20 , 20 ) , SED( 20 , 20 ) 

SU>U=0.0 

SUM2-0.0 

WIl«(N-M)/2fl 

NUTsN-lJNl+1 

DO  10  IsNNl.NNT 

DO  10  JtJINl.HUT 

10  SUI-ilsSUTIl+RL( l)*RL( J)*COV(I,  J ) 

DO  12  IsNNl.KNT 

12  SW2=RL(I)  *SX?(N,  l)+SUK2 
FKA«SUM2-SK 
FRA2sFNA*FKA 
SED(N,l!)=SUMi+FNA2 
RETURN 
END 


SUBROUTINE  PRINT( N , M. SEL , COFL, SSS , COFS , A , 13 ,  D) 
DIMENSION  SEl( 20.20)  ,COFL(20) ,SES(20,20)  ,C0FS(20) 

Nr=S AMPLE  SIZE  ' 

UPSIZE  AFTER  CENSORING 
SEL=KSE  LOCATION 

COFLsCOEFFICIENTS  FOR  LOCATION  ESTIMATE 
SESsMSE  SCALE 

COFSsCOEFFICIENTS  FOR  SCALE  ESTIMATE 
AsCONSTANT  FROM  LOCATION  ESTIMATE 
DsCOKSTANT  FROM  SCALE  ESTIMATE 
IF(I3,GT.5)  GO  TO  5 
I PAGE=? 1 
PRINT  20 

PRINT  21 , N , K. SSL( N , M) . 13 » COFL( 3) » 5ES( N , K) » 13* COFS( 3) 
PRINT  23, A, D 
X=38 
RETURN 
5  RM=5 
RK=K 

RLINSsRKf2 

IF(RLI!3E«GS,RX)G0  TO  10 
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25  AS ( I , J )=COV( 1+2 , J+2 )+EXP( N , 1+2 ) *EXP ( N , J+2 ) 

C  FILL  THE  BS  MATRIX 

DO  26  1=1, F4 
■V  \Z 6  BS(I)=EXP(N,I+2) 

C  COMPUTE  THE  COEFFICIENTS 

NTs=H 

CALL  MTXEL( NT , AS , BS , XS ) 

DO  27  1=1, NT 

27  RKC(I+2)=XS(I) 

CALL  SELD(  N ,  M,  1 , 0 ,  RXC ,  EXP ,  CO V ,  SES  ) 

SUMDbO.O 
DO  UO  ID=3»NP 
UO  SUMD=S UHDfRKC (ID) 

RKDsSUMD 

CALL  PRINT  (  N ,  M,  SED,  RL ,  SES ,  RKC ,  RK ,  3 , RKD) 

IF(N,EQ,5)  GO  TO  100 
IF(K.EQ.o)  GO  TO  100 
XF(N.EQ.7)  GO  TO  100 

C  CENSOR  AND  COMPUTE  THE  COEFFICIENTS 
GO  TO  29 

0  IF  IT  IS  DESIRED  TO  CENSOR  TO  SIZE  M=l,  REMOVE  PRECBEDING  CARD 
IF(iN/2)*2-N)28,29,28 

28  NC=(N/2)-2 
GO  TO  30 

29  HC=i:/2-3 

30  DO  50  IC2N=1,!!C 
NPCkK+1 
UP=K+2 

K=M-2 

DO  31  JD=3.NPC 

31  A(2, JD)=A(2, JD+1) 

DO  32  ID^,!JP 

DO  32  JD=3,NPC 

32  A(ID,  JD)=A(ID,  JDf  1 ) 

DO  33  ID=3.NPC 

DO  33  JD=2,KPC 

33  A(ID,  JD)=A(IDf  1 ,  JD) 

NP=Mf2 

CALL  MTXSL(I5P,A,B,X) 

DO  35  1=3, NP 

35  RL(I+ICEN)=X(I) 

RK=X(2) 

CALL  SELD(  K,  M,  RX,RL,  EXP,  COV,  SSD) 

C  COMPUTE  THE  SCALE  COEFFICIENTS  FOR  CENSORED  SAMPLE 
NTSsMfl 
NT=Mf2 

DO  36  IS =2, NT 
DO  3^  JS=1,IJT3 

36  AS(I$, JS)=AS(I3, JS+1) 

DO  37  IS=1, NTS 

DO  37  JS=1 , NTS 
BS(IS)=BS(ISa1) 

37  AS(IS, JS)=AS(IS+1, JS) 

CALL  l,TXEL(M,ASfL3,XS) 
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8  PRINT  21fK,M,SEI.(NlM),I3fCOFL(I3),SES(N,M),I3fCOFS(I3) 
M4ssI3fi 

if(m,eq,i)go  to  9 

PRINT  22,(1,  COFL(l ) ,  I , COFS ( I) , L-M4, M2) 

9  PRINT  23, A, D 
XsK-K-2 

GO  TO  16 

10  ISKIPsK+3 

DO  11  IS«1,ISKIP 

11  PRINT  24 
PRINT  25,IPAGE 
IPAGEsJPAGSfl 
PRINT  20 

K=4l 

15  GO  TO  8 

16  IF(N.£Q.20)G0  TO  17 
RETURN 

17  IF(M.EQ.1)G0  TO  18 
RETURN 

18  lSKIPssK+3 

DO  19  IS=1 ,ISKIP 

19  PRINT  24 
PRINT  25.IPAGE 
RETURN 


20 .  JSSSK1!?  •  i3HGl,r/.‘vlT!!/72-3>  / // , 41X, 8KTA3L3  II,//,19X,51HC03?F 

FOR  BEST  CONDITIONAL  ESTIMATION  OK  THE , / , 1?X , 5 -HLOCATION  A 
2®““*  PALTERS  OK  TriE  CAUCNY  DISTRIBUTION , / , 2oX , 3?N(aI TN  ADDI 

41 6?1,fi'wS1^u  .//» 15X ,  59( *.’*•),/ ,  23a  ,  14H**  LOCATION  *+ 

‘  V'^^lfi  3X,  1HM,  5X,  3NM3S, 5X,  1HI.4X,  5MC0BF, ,  10 

5X,3.tU'2>3,5X,lHI,4X,5HC0SF.,/,l5X#59(*.*)) 

21  F0RfIATClH0,l4xlI2,l4,Fl0.5,l4,Fll*6,F13.5tl4,?11.6) 

22  FORMAT (  33X ,  12 ,  FU ,  6 , 1  % ,  12 ,  Fl  1 . 6) 

23  FORI-IAT(34x,lHAtFll , 6, 16X,  1HD,F1 1 . 6) 

24  F0IUIAT(1H  ) 

25  FORMAT (42X, 13) 

END 


wppj'j*  j  wj'jl  4K1  pjwjJTJfiTWTfW 
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